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Abstract 

We summarize the present status of lattice gauge theory computations of the leptonic 

decay constants fo and Jb- The various sources of systematic errors are explained in 

a manner easily understood by the non-expert. The results obtained by the different 

. , groups are then systematically compared. As a result, we derive estimates for fu and 

[~^ I Jb in the quenched approximation through an appropriate combination of the data 

^ ' available from the different groups. Since we account for a possible lattice spacing 

dependence, the final errors are quite large. However, it is now well known how these 

O . uncertainties can be reduced. 

2^ ! For the decay constant of heavy-light pseudoscalar mesons with masses of 1-2 GeV, 

-)— » 



an interesting comparison of a full QCD result with the corresponding simulation in 
the quenched approximation can be done. Effects of sea quarks of mass rris are below 



p ■ the statistical accuracy of these simulations. 

^ ■ Related quantities, like ^-parameters, the spectrum of beauty-hadrons and the break- 

ing of the QCD string are discussed briefly. 
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1. Introduction 



During the past decades, experiments in high energy physics have resulted in quite 
an indepth understanding of the interactions between quarks and leptons. It is math- 
ematically formulated as the Standard Model of particle physics. The fundamental 
parameters that characterize the overall strengths of the electromagnetic, weak and 
strong interactions in the Standard Model have been determined precisely in experi- 
ments. Also, the masses of the matter fields - apart from the postulated Higgs-boson 
and top-quark - are known. However, the detailed structure of the weak interactions 
of quarks is not determined well. Let us discuss this in more detail. 

In the Standard Model, quark fields q couple to the charged weak interaction vector 
bosons iy+ through a term 

^ Y. ^u 7m — Y^ VuD Qd W^ . 

U=u,c,t D=d,s,b 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix|T| Vud originates from the trans- 
formation from weak interaction eigenstates to mass eigenstates. Assuming that there 
are only three generations, it is unitary and can be written in terms of four observable 
parameters. 

Given the strength of the weak interactions through fi - decays, in principle one needs 
to determine the rate of decays of quarks such as d ^ u to obtain the matrix elements 
of V. However, quarks do not exist as free particles and we can only observe the 
decay of hadrons. In general, the full knowledge of the hadron wave function is needed 
in order to relate a measured decay rate of a hadron to the parameter Vud in the 
Lagrangian. 

Fortunately, the neutron and the proton are related by isospin symmetry which is well 
tested experimentally. Hence, the vector part of the decay amplitude of the ordinary 
P decay {n —>■ p) is given by a Clebsch-Gordan coefficient times the desired parameter. 
So \Vud\ can be determined with almost no theoretical input. 

Analogously, \Vus\ can be extracted from K ^ it decays starting from an approximate 
S\J{3)fiavor symmetry. Here, it is already important to include symmetry breaking 
corrections, which can be done reliably using chiral perturbation theory [Q. 

With these two matrix elements determined from experiment|0] ( \Vud\ = 0.9744(10), 
\Vus\ = 0.2205(18) ) and imposing the constraints of unitarity, Wolfenstein observed 
that mixing appears hierarchically ffl. The hierarchy is parametrized by a small pa- 



rameter A = IK, 



us] 



Vud Vus Vub\ I 1-AV2 A AX\p-i7]) 

V=\ V,d Vcs VcA = \ -A 1 - AV2 - lA^X^ A\^ 

Vtd Vts Vtb J \A\^{l-p-tr]) -AX^ 1 

The other parameters A, p and rj are of order 0{X^). 

The parameter A has been determined from semileptonic decays B -^ D*lv. As 
theoretical input, one needs some information on the formfactors of these transitions. 
This has previously been extracted using model wave functions for the B and the 
D-mesons 0. It has been discovered, however, that the theoretical description of 
these formfactors is considerably simplified in the limit when the masses of both the 
h and the c quark are large compared to typical hadronic scales. In this limit, there 
appears an approximate SU(2)//a„orXSU(2)5pj„ symmetry relating the 6-quark and 
the c-quark. Consequently, the number of formfactors describing these transitions is 
reduced to one, the Isgur-Wise function]^. The effective theory that starts from this 
symmetry and tries to include the corrections of order 0{mN/'mh), with rrih the heavy 
quark mass and rriN the nucleon mass representing a typical QCD scale, as well as the 
QCD radiative corrections, is called Heavy Quark Effective Theory (HQET). Apart 
from radiative corrections0, the normalization of the Isgur-Wise function at the point 
v.v' = 1 [v and v' are the 4-velocities of the initial and final state meson respectively) 
is given by the symmetry up to corrections of order 0{{mj^ /rric)'^) 0. Assuming that 
the latter terms are small, one extrapolates the experimental data to that kinematical 



point and one obtains |llO|, p!T| A = 0.90(12). 



Again, a symmetry has helped to circumvent the full solution of the bound state prob- 
lem of QCD. It is important to note, however, that in the case of the HQET, we have 
to date a far more limited understanding of the size of symmetry breaking terms than 
in the case of S\J{3) flavor- In the latter case, there is a large amount of experimental 
information on symmetry breaking, which helped to develop a comprehensive theo- 
retical treatment of these effects. Therefore, it is of interest to quantify or bound the 
0{{'mp/mc)'^) corrections in the above analysis. A promising approach to this prob- 
lem are lattice gauge theory calculations of the semileptonic form factors[0, |13|. In 
addition, QCD-sumrules allow to estimate the size of the 0{{mp/mc)'^) corrections |14 



There are two additional unknowns in the CKM-matrix: p and rj. They are of particu- 
lar interest, since they are a (and probably the) source of CP-violation in the Standard 



Model|]l5[. Their values are mainly constrained by three experimental observations. 



i?-meson decays together with model calculations determine ||rT| y/p'^ + rj'^ = 0.36(9) 



Moreover, the CP-violation parameter |e| in the /r-system [p!5|] and the Bq-Bq mixing 
parameter x^ restrict p and rj. However, in order to extract p and rj from these 



measurements, one needs to know the mixing matrix elements 



3 

and 



i?B/|M| =< Bolid^^il - 75)&)(rf"7M(l - 75)&)|5o > 



where Jb and Jk are the leptonic decay constants of the B- and the i^-mesons and Bb, 



Bk parametrize the matrix elements relative to the vacuum insertion "approximation" |jl5| 



The dependence of the allowed domain in the (p, r/) plane on the theoretical uncertainty 
of Bk [Jk is known from experiment) is not very large. The allowed region depends 



sensitively on the value of the product -B_b/^, however |lq, |T^ |n 



Until fairly recently, phenomenological analysis of the allowed region in the (p, t]) plane 
used the "old prejudice" fsy/Bs = (110 — 160) MeV. Such values favor a negative 
p and lead to small predictions of CP-asymmetries in B decays, which one would 



like to measure in future experiments at -B-factories]^. First lattice estimates [|1^, ^ 
of these quantities that were done within the HQET approximation, indicated much 
larger values of Jb- Also, QCD sum rule calculations performed in the HQET limit 



have subsequently yielded such large values|2y, ^. Together with Bb — 1 which is 
supported by a number of theoretical investigations, this indicates that positive values 
of p are favored by the experimental constraints. If this is the case, much larger signals 
for CP-violation in B decays are expected - an interesting perspective for i?-factories. 

We emphasize, that the value oi BbJb is the central theoretical uncertainty on our way 
from the experimental data to the determination of the CP- violating part of the CKM- 
matrix. It is therefore important to determine BbJb with a good precision in order 
to narrow down the allowed region of {pyt]). To this end, precise lattice calculations 
of these quantities are necessary and one has to take account of the full spectrum of 
the statistical and systematic uncertainties. Besides this number (and similar other 
ones that are of interest), lattice gauge theory calculations also offer the possibility 
to quantify the corrections to the HQET limit. Starting from the simpler case of the 
D-meson, we review the present status of these efforts and their future perspective 
in sections ^, ^ and [6]. A discussion of some general features of lattice gauge theory 
calculations - as given in section - is needed to understand the different sources of 
systematic errors in the computations. Other sections cover related topics. 



2. Lattice QCD 



Lattice QCD[^, |2^ is a regularization of the fundamental theory of strong interactions 
that allows nonperturbative calculations. In principle, it involves no approximations. 
The latter arise because our numerical possibilities of evaluating the path integral by a 
Monte Carlo process are limited. This means that computations have to be performed 
at parameters that are different from their values in nature. Therefore, extrapolations 
have to be performed e.g. in the lattice spacing, quark masses and the space-time - 
volume V. Furthermore, we are limited by the statistical errors of the Monte Carlo 
calculation. 

In order to be able to discuss the various systematic errors in current calculations, 
we need to introduce some definitions and summarize the crucial points in a practical 
lattice QCD computation. This section primarily addresses readers that are not very 
familiar with lattice QCD. The specialist may, however, be interested in the discussion 
about universality in the quenched approximation, which is given in section B^B. 



2.1 Path Integral Representation of Greens Functions 

The starting point of lattice QCD is to define the field theory on a space-time lattice 
with spacing a. This is one way of regularizing the ultraviolet divergences that are 
intrinsic to a quantum field theory: the lattice spacing serves as a cutoff for the high 
frequency modes of the fields. As in the case of any cutoff, the field theory is finally 
defined in the limit a~^ —>■ oo. 

The basic variables of lattice QCD are the dimensionless quark fields qf{x) (/ labels 
the different flavors of quarks and we omit the color- and Dirac-indices) and the gluon 
field variables U^{x) which are SU(3) matrices. The coordinate x resides on the 4- 
dimensional space-time lattice: x = {xq,Xi,X2,X3), x^/a e N. The gluon variables 
U^{x) are the gauge connection of nearest neighbor points x and x + afi. 

The dynamics of the theory (i.e. bound state masses and properties, transition proba- 
bilities) is contained in the euclidean Greens functions or correlation functions. They 
are defined through the Feynman path integral 

<0[g,g,f/]> = Z-^ j D[q] j Y)[q] j Y)[U]eM-S[qA.U]) 0[g,g,f/] 

Z = JDlq] JDlq] J D[U]expi-S[q,q,U]). (2.1) 

Here, / D[g] / D[g] stand for multiple Grassmann integrals over each component of the 
quark variables and / D [U] represents an integral over the gluon variables on each link 



with the SU(3) Haar measure. 0[g, g, U] is a functional of any number of variables of 
the theory and S[q, q, U] is the euclidean action 

S[q,q,U] = SG[U] + Sp[q,q,U] (2.2) 

with a gluonic part 

Sg[U] = /5E E {1 - i^e tr P.^x)} ; /? = -| (2.3) 

and a fermionic part 

SF[q,q,U] = Y.^{x)Vq{x)- (2.4) 






ycr,, 



Here, Pf^^^ix) is the product of gauge parallel transporters f/^(x) around an elementary 
plaquette 

PiJL,v{x) = U^{x)Uy{x+ajl)U^^{x+ajl+ai))U^^{x+ai') ; f/_^(x) = Uj^{x—ajl) , (2.5) 

the covariant finite differences are 

V;.g(x) = U^{x)q{x + afi) - q{x) , V; = -V_^ (2.6) 

and the euclidean gamma matrices fulfill 

{-f^^,-fu} = 2Sf,,u , 7^ = 7/., (r^^,y = -[-fi,,lu]- (2.7) 

A lattice approximation of the field strength is given by 

^/.^(^) = g [ Pf^Ax) + P-,^,-,.ix) + Pu,-f,ix) + P-uAx) 

- P,,^(x) - P„.,-^(a;) - P_^,.(x) - P^,-.(x) ] . (2.8) 

The lattice action S depends on the bare parameters of the theory, namely g^, the gauge 
coupling, and the hopping parameters Kf that are collected in the bare quark mass 



matrix Mq = diag(l/(2i^i), ... , l/(2i^„ ))— 4. The last term in the fermion matrix [p^ 
V, which is accompanied by the coefficient function csw{9o) "^i^^ be discussed at the 
end of section W^ 



Eq. ( p.l|) looks identical to a statistical mechanics thermal average. It can be evaluated 
by Monte Carlo importance sampling, once the space time grid is restricted to a finite 
number of pointsp3[. Such a Monte Carlo evaluation is the main nonperturbative 



method employed in lattice gauge theory calculations. There is an essential difference 
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to a simulation of a spin model, however, namely the fermionic nature of the quark 
fields. Since the fermionic action is a quadratic form, the Grassmann integrals can be 
performed analytically, e.g. 

/D[g] |D[g] g^:(x')g^(x) = p-},„,,^,,^,, ,,,,[[/] deiV[U] eMScP]) , (2.9) 

with a, a' Dirac indices and c, c' color indices. In general, one obtains 

< 0[q,q,U] >= Z-^ fD[U] expi-Sap]) detV[U] Fo[U] (2.10) 



with a functional -Fb[?7] that is easily derived using eq. ( p.9|) and the property that 
the Grassmann fields are anticommuting. It is eq.( |2.1CI| ) which is evaluated by Monte 



Carlo important sampling, generating "configurations" U with the probability P[U] oc 
exp{—Sa[U]) detV[U] and then averaging -Fb[t/] over these configurations. 



2.2 A Pseudoscalar Correlation Function 

We discuss here the most simple (and for the purpose of this review almost suffi- 
cient) observable that is well calculable in lattice QCD. We choose a lattice with a 
finite number of points L/a in the three space directions and Lt/a points in euclidean 
time. The gauge fields are periodic functions: Ufj^{x) = U^{x + kL), k = 1,2,3; 
U^{x) = U^{x + OLt). The quark fields are taken antiperiodic in time and periodic (or 
antiperiodic) in space. We may take the axial vector current 

Af/{x) = Za qfixh,l,qf'ix) (2.11) 

as a convenient interpolating field for pseudoscalar mesonsQ 



<0\Ao{0)\P>=FpJMp/2 . (2.12) 



The finite renormalization Z^ of the axial vector current will be discussed in section |27 
The (space-) momentum zero correlation function of Aq has a spectral decomposition 

J2 < A{x)Al{0) >= E I < 0\Ao\P,n> p{exp[-E„xo]+exp[E„(xo-Li)]} , (2.13) 
X "-^ 

where \P,n > denote the excited states in this channel and \P, 1 >= \P >. 

With the Wilson action (i.e. cswiOo) — 0)) ^^^ transfer matrix, which describes the 
propagation in the euclidean time, is positive and hermitian |p5|| . In this case, eq. (|2.13|) 



is therefore exact up to irrelevant terms of order O {exp[— EgapLt]) , where Egap is the 

^ |P > denotes the pseudoscalar state with momentuni zero, the appropriate flavor quantum 
numbers and with a normalization < PIP >= 1. 



lowest state with the quantum numbers of the vacuum. When cswiOo) 7^ 0, one 
expects eq. ( p.l3|) to hold to a good approximation at moderately large Xq. 



At sufficiently large values of the euclidean time xq, the spectral representation shows 
that the correlation function eq. ( p.l3| ) is dominated by the lowest state and one may 
extract both the mass Mp = Ei and the leptonic decay constant Fp from a Monte 
Carlo estimate of < Aq{x)Aq{Q) >. 



2.3 The Continuum Limit 

In this section, we discuss how the continuum limit a — >^ is taken in QCD. 

Lattice gauge theories contain only dimensionless variables and fields. All quantities 
are calculated in units of the lattice spacing a. The latter is not a parameter of the 
calculation but is determined aposteriori once a dimensionful observable is fixed to its 
experimental value. 

In order to shorten the discussion, consider QCD with only mass-degenerate quarks, 
i.e. Kf = K for all /. Eq. ( |2.13| ) then determines the dimensionless functions]^ 

Fp{9lK) (2.14) 

and 

Mp{glK). (2.15) 

In addition, for example, the mass of the nucleon in lattice units M]\f{gQ,K) may be 
calculated from a nucleon correlation function. 

Since the quark masses are free parameters in the Standard Model, we have to set a 
nonperturbative renormalization condition to fix the hopping parameter K in terms 
of a renormalized, experimentally observable quantity. A convenient choice is 

Mp{gl K)/Fp{gl K) = mJU , (2.16) 

where /,r and m.,^ refer to the experimentally measured leptonic decay constant and 
mass of the pion. Eq. ( |2.16| ) defines a renormalization group trajectory in the [g^, K) 



plane. Along this trajectory, the ratio 



Fp{9lK) 



Rmn,Fp{.Fp) = ^ ,^2\^\ \Mp{gl,K)/Fp(glK)=m^/U {'^■^^) 



^ We assume that at each point {g^ , K) , the space extent L of the lattice has been chosen large 
enough such that finite size effects are negligible compared to the statistical accuracy of the calcula- 
tion. Therefore, the dependence on the third dimensionless parameter (.Fp L/a) is not written. 
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is a function of Fp only and reaches its continuum limit at Fp = since Fp is the 
physical decay constant times the lattice spacing. From asymptotic freedom one infers 
that the continuum limit is reached at the point g^ = 0, K = 1/8. 

Equivalently, on the trajectory eq. (|2.16|) , we may consider the lattice spacing to be 
determined by 

afM) = Fpi9lK)/U (2.18) 

and the nucleon mass in physical units is given by rriM = M^igo, K)/af^{gl). 

Corrections to the continuum limit constitute a major source of uncertainty for prac- 
tical lattice gauge theory determinations of quantities like fs- These corrections have 
to be addressed by Monte Carlo simulations done at different points on the renor- 
malization group trajectory. Nevertheless, it is instructive to consider perturbation 
theory and learn something about the structure of the finite lattice spacing correc- 
tions. Symanzik has discussed the cutoff dependence of Feynman diagramsp6|. For 



the action eq. (|2.4|), such a discussion ||2^ suggests that the leading dependence of phys- 
ical quantities such as Rmi^,Fp{,Fp) on the lattice spacing is linear (up to numerically 
unimportant logarithmic modifications): 

Rm^,fAFp) = ^ + PiFp + 0{Fl) ■ Fp = Uaf^ . (2.19) 

How small does the lattice spacing have to be such that the finite a corrections are 
small? Again this can only be answered by a nonperturbative calculation of the func- 
tion RMf^,Fp{Fp) and for each such a function separately. It is however worth noting, 
that one may expect that the relevant scales are the sizes of the wave functions of 
the hadrons considered (for S-wave states), not their Compton wave length, since we 
mainly need a good approximation of the continuum wave function by the lattice one. 
We may therefore hope that eq. ( |2.19| ) can be applied for a « 0.5/m. 



Improving the Approach to the Continuum Limit 

In a formal expansion in powers of the lattice spacing, the fermion action eq. 
breaks up into three pieces. The first two, the naively discretized Dirac operator and 
the Wilson term, are necessary to obtain the correct continuum limit [^. The Wilson 
term, ^S/^ V^V*g has a classical continuum limit —aqAq + 0{a^), where A is the 
covariant Laplace operator. This dimension 5, irrelevant, operator is responsible for 
the lattice spacing corrections of order 0{a) as e.g. in eq. ( |2.19| ). It is expected that 
for on-shell matrix elements these terms can be cancelled by adding one dimension 5 



operator to the action|28|. Most conveniently, one chooses the third term in the action 
with coefficient function cswigo)-, since it does not involve further derivative terms in 
the quark fields p^. 



Symanzik's idea was to calculate coefficients like cswigo) ^^ perturbation theory lead- 
ing to a perturbatively improved action. Already the tree level expression cswido) = 1 
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leads to pi = O^g^). Since i) the arguments that lead us to expect that the 0{a) lat- 
tice artifacts can be cancelled in this way are due to perturbation theory and ii) 
the coefficient function cswiGo) is calculated perturbatively, improvement needs to be 
demonstrated through Monte Carlo results. At present, the evidence for improvement 
with the action of Sheikholeslami and Wohlert is somewhat indirect but nevertheless 
suggests that 0{a) lattice artifacts are significantly reduced compared to the Wilson 
action|2^. In this context, it is furthermore of interest that in the case of the pure 
gauge theory with Dirichlet boundary conditions, an improvement has been estab- 
lished unambiguously for a one loop Symmanzik improved action||5^. Thus, there is 



an example which shows that the perturbative improvement program of Symmanzik 
works at those values of the gauge coupling that are needed in realistic simulations. 

Below, we will also discuss recent results for heavy-light (HL) mesons that were ob- 
tained with the perturbatively improved action. We will use the synonyms "improved 
action" and "SW-action" for the action eq. ( p^ ) with cswigo) = 1, while we use the 
name 'Wilson action" when csw{go) = 0. 



2.4 Renormalization 

In the previous section, we have ignored the problem of the determination of the 
strength of the axial vector current Za- At first sight, it requires a nontrivial renor- 
malization since it is a composite operator. However, at the formal level, the flavor 
offdiagonal axial vector current A^'^ (x) is the current of chiral symmetry, which is 
broken only softly: 

d^Af/{x) = {rrif + mj,) P^'^'(a;); P^^^ (x) = qf{x)^^qf>{x) . (2.20) 

This relation (more precisely the chiral ward identities) would normally insure that 
A^j, does not get renormalized. The quantum field theory needs to be defined with a 
regulator, however. As the lattice regularization (eq. ( p.4|) ) breaks chiral symmetry by 
terms of order 0{a), the formal argument does not apply and A^ acquires a nontrivial 
renormalization |]3T1, |32|.P| 

In the following, we mainly discuss the case of the Wilson action, which is still the 
action for which most results exist. At the very end of this section, we mention the 
case csw = 1- 

We write the relation between the bare lattice current and the renormalized current 
Ajj^f (x) in the most general way: 

4'^'(x) = ZA{gl Kf, Kj,) qfixh,^,qf,{x) . (2.21) 



Only in the case of the vector current, a conserved current can be defined for most actions. 
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Here, Za needs to be defined through an appropriate normalization condition, such as 
eq. (|2.27|) which will be discussed below. We may rewrite eq. (p.21|) in a form that was 
suggested by Lepage, Kronfeld and Mackenzie p3|, |3^ , |35| P| 



A^f{x) = ZAigl Kf, Kf,) J^^ - 3u J^^ - 3u qf{x)^,^,qj>{x) , (2.22) 

where the mean field value u of the gauge field is defined through a gauge invariant 
quantity, e.g.[p^ 

U = l/{SK,r^t) . (2.23) 

The critical value of the hopping parameter Kcrit is to be determined nonperturbatively 
from the Monte Carlo calculations through Mp{gQ, Kcrit) = 0. 

If the renormalization constants Z^ and Z^ are defined through a nonperturbative 
normalization condition at a finite value of the lattice spacing, eq's. ( p.21| , ^^22D are 
completely equivalent. In most applications, the normalization of the currents is taken 
from perturbation theory and after dropping the 0{a) terms. In this approximation, 
there is a significant difference between the currents ( ^.21| , pT^^ . Numerically, the 
difference is dominated by a factor exp[(m/ + m//)a/2]; anif = log[l + {2Kf)~^ — 
[2Kcrit)~^]- Although this factor becomes irrelevant in the continuum limit, it is 
sizeable in present HL calculations. 

In 1-loop order perturbation theory one obtains [^, ^, ^ 

ZAigo,Kf,Kf') = 1- 0.133373^^ + 0(^^) + 0(am/) + 0(am//) (2.24) 
ZAiglKf.Kf,) = I - 0. 02A80 g^ + O (g^) + {am f)+0 {am f) (2.25) 

~g^ = gl/P \ P=\< trP,A^) > . (2.26) 

In the above equations, we did not need to specify the normalization condition for the 
current, because the Z-factors depend on these conditions only when 0(ar7ij)-terms 
are included. Instead of using the bare coupling as expansion parameter, we have 
written the perturbative expansion of Za, Za in terms of the effective coupling g"^, 
introduced by G. Parisi|3^. Through several examples, Lepage and Mackenzie [p7| . 



demonstrated that the series in g^ (at momentum scales of the lattice cutoff vr/a) 

^ The argument that is given in ref. |33t| for this normahzation is based on the hmit auif >> 1 
and is therefore only indicative in the present context. On the other hand, Kronfeld p4| s tarts from 
the canonical normalization of the quark fields in the transfer matrix representation||25| of Green 



functions like eq. (2.13). Subsequently, he uses a mean field approximation and implicitly assumes 
that the quark fields enter with zero space-momentum. One is therefore well advised to consider this 
normalization as an alternative and not as the solution to the problem of handling heavy quarks with 
masses of the order of the inverse lattice spacing. 
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converges more rapidly than the one in the bare couphngQ . In the known examples, 
the correction terms to the 1-loop order are compatible with having coefficients of the 
order one or less, when the series is written in terms of a = g'^ /{A-k). In addition, 
Lepage and Mackenzie present examples where 1-loop perturbation theory gives quite 
accurate results when compared to full nonperturbative computations of the quantity 
considered; see also ref. [^. 



It is argued in ref. [Q that eq. ( p.22| ) is the definition with the smaller lattice artifacts 
when it is used together with the perturbative value of the renormalization constant 
(as given above). In subsequent works, this normalization is quoted as the "correct" 
one and eq. ( p.22| ) as "false" . We stress that this is not the case as a general statement 
and give a counter example in the appendix. The example, as well as the lattice 
spacing dependence of /o, which we investigate later, mainly teaches us that the size 
of the lattice artifacts has to be determined by precise simulations with varying lattice 
spacings for each case separately (A one-loop calculation including the 0(m/a)-terms 
could shed some light onto this question as well). 

The 0{g'^) corrections to the one loop estimates of the renormalization constants are 
a systematic error. In principle, the axial current may be renormalized nonpertur- 
batively. Here we only outline the idea and refer the reader to ref. [^ for a more 
detailed discussion. The idea is to require current algebra relations to be valid for the 



renormalized currents. For instance, one may require [p2 



^ X 

= -{Vj'-f'{y)Vj''f{z)) + {Af/{y)Al''-f{z)) (2.27) 



for the renormalized currents. Taking u to be i) spacelike and ii) timelike, Eq. (p.27|) 



constitutes two rather independent linear equations for the current renormalization 
constants. As these equations are nonlinear, the renormalization constants can be 
computed - without the use of perturbation theory. 

It is important to note the following points about this equation^. The axial vector 
current at point y on the left hand side is changed into the vector current on the right 
hand side by the contact with the pseudoscalar density and analogously for V{z). This 
means that the relevant contributions in the integrated three-point function on the 

^ Lepage and Mackenzie actually go further pa and change from g^ to a physical coupling defined 
from the interquark potential. In addition, they give a prescription how to optimize the scale at which 
this coupling should be evaluated when one inserts it into a 1-loop perturbative result. Although 
this does indeed improve the agreement of l~loop perturbation theory with nonperturbative Monte 
Carlo results for the examples that are considered in p5| , it is important to be aware that the errors 
remain of order 0{a^). When perturbation theory is used in this review, we use the expansion in 
terms of a = g^/(47r) and try to indicate the expected size of the 0{a^) terms. 

^ We would like to thank M. Liischer and S. Sint for pointing out this problem. 
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left hand side originate from short distances \x — y\ and |x — z|. At distances of a 
few lattice spacings, lattice correlation functions do exhibit a strong dependence on 
the lattice spacing. This will translate into potentially large lattice artifacts of the 
renormalization constants. A sufficiently small lattice spacing is therefore required 
to be sure that the 0(a)-terms injected through eq. ( |2.27| ) are numerically small. At 



current values of the lattice spacings, it is not evident that this requirement is fulfilled. 

One can avoid this problem at the price of giving up one prediction: In ratios of 
decay constants, the renormalization cancels and only unavoidable 0{a) terms remain. 
Unfortunately, this option increases the statistical errors in the decay constants fo and 
fs and consequently the remaining 0{a) terms are more difficult to control. 



We point out that eq. (|2.22|) is only an ansatz. In order to systematically reduce the 



0{amf) corrections one has to both improve the action through csw 7^ and the 



operators. With csw = 1, one can achieve 0{amf) -^ Ojaamf)'^^ P9l 



In summary, the renormalization of the axial vector current is a nontrivial problem 
in practice and introduces an uncertainty in the determination of pseudo scalar decay 
constants. The same holds true for the normalization of four fermion operators such as 
the one responsible for BB mixing. Furthermore, when working in the HQET in order 
to evaluate the limiting behavior of the decay constant as one of the quark masses 
becomes very large, the current develops an anomalous dimension and - at present - 
no practicable nonperturbative method to renormalize the current is known. 



2.5 Reaching Asymptotics 

In a Monte Carlo simulation one evaluates correlation functions such as eq. ( |2.13D at 
finite values of the euclidean time xq. In order to be able to extract hadron masses and 
decay constants, one needs to insure that the terms with n > 1 are negligible in the 
spectral representation eq. ( |2.13| ). For a quantitative investigation of this problem, we 
consider an example that is relevant for the physics discussed in this review, namely 
the correlation function of a heavy-light (HL) axial vector current, composed of one 
heavy quark (with a mass of the order of the charm quark mass) and one light quark. 

It is useful to look at the logarithmic slope of the correlation function. This slope, 
called "effective mass", approaches a constant when the above mentioned corrections 
are negligible. In fig. |l] (filled symbols), one observes that the effective mass shows a 
slow variation with the euclidean time still at a distance of xq = 15a ~ 5GeV~^. At 
larger xq, a plateau seems to develop in the effective mass. The region of the plateau 
in fig. |1| is rather short, however. 
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Figure 1: The logarithmic slope of the pseudoscalar correlation functions C'°'^''°'^(xo) (filled symbols) 
and C'^''^{xo) (open symbols) |Q. G denotes a gaussian wave function with r.m.s. radius of about 
0.3 fm. The cutoff is a~^ = 3.17GeV (/3 = 6.26) and the simulation was done in the quenched 
approximation. The light quark mass was chosen about twice the strange quark mass (K=0.1492), 
while the heavy quark mass is roughly the charm mass (K=0.1350) in the lower part and somewhat 
higher in the upper part (K=0.12). The dashed curves illustrate the effect of the finite time extent 
of the lattice according to eq. (^.3l|), assuming that only one state contributes in this range of a;o- 
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One may improve the situation by choosing better interpolating fields for the hadrons. 
We briefly describe the approach that has been taken so far for the case of meson 
flelds. The general idea is to "smear" the local quark flelds with trial wave functions 

q\x) = Y. ^\x, y; U{yo)) 6^,y, q{y) . (2.28) 

y 

An interpolating fleld for a meson is then given by 

M'{x) = q{x)^o^^q\x) , (2.29) 

where the index J denotes the trial wave function and we have omitted the flavor 
indices^. $*^(x, y; U{yo)) depends either in a gauge covariant way|]4l|, ^, ^, ^ on the 



link variables denoted by U{yo), or the construction is done in Coulomb gaugeH^. A 
meson correlation function is constructed by 

C''^{xo) = E < -^'(^) [-^-^(0)]^ > • (2.30) 

X 

It has the spectral representation 

^''■"(^o) = E < 0|^'(0)|P,n X P,n\[M'{0)]^0 > {exp[-EnXo]+exp[En{xo-Lt)]} 



n>l 



(2.31) 

In fig. we have included the logarithmic slope for such a smeared correlation function [44 
It reaches a plateau at values of xq ~ 5a - much earlier than the correlation function 
of the local axial vector current. 

A clean test for ground state dominance is provided by the ratio 

R (^o) - [^/,.(^^)]2 ' (2-32) 

which is constructed such that it becomes one when all correlation functions are dom- 



inated by the ground state|4^, ^. This is a very stringent test since contrary to the 



search for a plateau in the local mass, the value of the plateau is known here too. 

Since most computations of Jd have been carried out implementing the local axial 
vector current, it is of interest to use this ratio to further investigate in how far the 
correlation function is dominated by the ground state at accessible values of xq- So 
we display R^'''°'^{xo) in fig. ^ where G refers to the gaussian smearing already used 
in fig. U and loc denotes the local axial vector current. R'^'^°'^{xq) is only consistent 
with one for the largest values of xq- Therefore, we consider the region of xq, where 



^$ is the wave function in terms of the relative coordinate. In addition to the form eq. (2.29), the 
case where both quark fields are smeared has been considered. 
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Figure 2: The ratio i?'^''°^(xo) is shown for a Ught quark mass of about twice the strange quark mass 
and a heavy quark mass of roughly the charm mass. Open circles are for a cutoff of a^^ = 3.17GeV 
(/3 = 6.26) and full points are for q ~^ — 2.32GeV (/3 = 6.00). The lattice spacing and thus the scale 
on Xq are defined through eq. ( ^.ij ). 



all correlation functions are dominated by the lowest two states. There, R ' (xq) is 
described by 



[1 + ej exp(-xoA)] [1 + e^ exp(-xoA)] 

l + e/ejexp(-xoA) 
<0|A^^(0)|P,2>/<0|A<^(0)|P,1 > , A = E2-Ei 



(2.33) 



Fig. D suggests that the smearing wave function was chosen well resulting in ec << Uoc- 
We may therefore simplify further: i?*^''°'^(xo) ~ [l + ef^^exp{—XoA)]. A fit to this form 
is shown as dashed curve in fig. 0. It describes the ratio well with a gap A ~ 600MeV. 
We have thus roughly quantified the dominant correction to the asymptotic behavior 
of the local axial vector correlation function. The size of this correction term shows 
us that some of the earlier lattice computations that investigated the dependence of 
the decay constant of HL mesons on the mass of the meson are subject to a significant 



systematic error since they used only local meson fields. One should not regard this 
investigation as a quantitative determination of the gap in this channel. 

If one uses C^°'^'^°'^ and does not isolate the ground state properly, one obtains an over- 
estimate of the decay constant. Repeating the above analysis at different values of the 
heavy quark mass shows that ef^^ increases when the heavy quark mass increases, while 



A stays roughly constant ^J]. Therefore, the effect becomes stronger with increasing 
quark mass, resulting in a misleading quark-mass dependence. Physically, this effect 
means that the ratio of the leptonic decay constant of the excited state to the ground 
state decay constant increases with the mass of the heavy quark. Unfortunately, at 
present, the precision is not high enough to determine ef^c and its mass dependence 
well. 

A further comment is in order concerning fig. ^ the figure contains results for two 
values of the cutoff. The ratio has a continuum limit only in the region of xq where 
it is approximated well by the truncated form R^'^°'^{xo) ~ [1 + efQ^exp(— a^oA)] that 
was used in the fit. At smaller xq where this is not true, nonuniversal matrix elements 
of the trial wave functions appear and there is no reason to expect an (approximate) 
independence of the lattice spacing. Therefore, the a-independence of the ratio that 
is observed for xq > 1.5GeV~^ is consistent with ec « eioc, that means a well chosen 



trial wave function in refs. [HSl K 



As a consequence of this investigation, it appears mandatory to use smeared meson 
fields in order to reach asymptotics and to perform a reliable computation of the decay 
constant^. In order to extract fp, the matrix element < 0\A4^{0)\P,1 > which is 
introduced through the trial wave function, must be cancelled. As seen from eq. ( |2.31|) 
this can be achieved by combining the results of an analysis of both C^'^xq) and 
C^'^"%xo). 



2.6 The Quenched Approximation 

As we mentioned in section (2.1), the first step in a Monte Carlo evaluation of the 
expectation values eq. ( p.l|) consists of generating lattice gauge fields U with proba- 
bility exp{—SG[U]) detV[U]. The only algorithm which generates this distribution 
exactly^ is the Hybrid Monte Carlo algorithm (HMC). It is numerically quite slow. 
In particular, the required CPU time for a given simulation scale^ approximately 
like {mga)"'^ with rj r^ 5 [Q. At present, this allows for exploratory simulations at 

^ Of course, the problem can also be solved by going to sufficiently large values of xq on a large 
lattice with large statistics in the Monte Carlo [Q. 

^Actually, the HMC algorithm[^ is valid for pairwise mass degenerate quarks only. 
^° This scaling law includes the scaling of the box length L that is necessary to avoid finite size 
effects when going to small quark masses. 

19 



unphysically heavy masses niq of the quarks. 

In perturbation theory, detX'[[/] generates quark loops. One may assume as a starting 
point, that for the full path integral at a given value of the cutoff, the effect of the 
quark loops is mainly to modify the renormalization of the coupling. Therefore, low 
energy observables may be approximately represented by a model with detX'[f/] -^ 
const. Also the success of constituent quark models suggests that neglecting quark 
loops in this way may be a sensible approximation to QCD. Clearly, we do expect 
significant differences between QCD and its quenched approximation, in quantities like 
flavor singlet densities in hadrons or the connection between the high energy running 
coupling and hadron masses. In this sense, it is better to think of the "quenched 
approximation" as a model. It does not constitute a systematic approximation to 
QCD and its accuracy depends strongly on the quantity that one considers. 



<3 

E 




mq/ms 

Figure 3: Extrapolation of hadron masses to the chiral limit. The extrapolation curves are used 
to locate the physical point in the quark mass, that means one point on the renormalization group 
trajectory. The ordinate mq/mg is the bare quark mass in units of the strange quark mass. mp{mn) 
denotes the p-meson mass at the physical point. The figure is taken from ref. [Q; error bars are 
smaller than the size of the symbols. 



Having emphasized this, we proceed nevertheless to summarize numerical evidence 
that the "quenched approximation" is surprisingly close to the real world. We consider 
the low lying hadron mass spectrum. The large effort that went into computing the 
hadron mass spectrum has been reviewed in references |^T], ^. A significant step on 
the way towards a reliable computation of the hadron mass spectrum has been taken 
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by the API5-group, using nonlocal trial wave functions (cf. section 2^) and quark 
masses around the strange quark mass[^ . 



This has recently been improved further in the GFll-project||5^ reaching a rather 
high statistical precision. The latter allowed for a more systematic study of the lattice 
spacing dependence and the finite volume effects. Using the assumption that up and 
down quark masses are degenerate, and the quenched approximation, the problem has 
two parameters only: K = Ku = Kd and g'^. To get to the physical masses in the 
continuum limit, one proceeds as discussed in section p?3| . An important difference to 
the idealization in section |2.3| is that one cannot compute the quark propagators at 
the physical values of K since this corresponds to the solution of an almost singular 
system of linear equations. Instead, all hadron masses are calculated at various values 
of the hopping parameter K and then extrapolated using forms suggested by chiral 
symmetry (for the pion mass) and simple mass perturbation theory. As an illustration, 
we show the extrapolations that were done in ref. |^^ in fig. |3|^^| . 



In ref. |Q the value of K at each Qq was fixed requiring Mv{K, gQ)/Mp{K, Qq) = 
mp/rriTr, where My is the vector meson mass in lattice units and rrip the experimental 
p-meson mass (the p-meson is a stable particle in the quenched approximation). 

The crucial point is now to calculate hadron mass ratios on the renormalization group 
trajectory and investigate their dependence on the lattice spacing measured as Mp = 
rUpa. Fig. H shows the example of ^f- {Mp) . A linear dependence on Mp (see section pl3|) 
is in agreement with the data and gives a value in the continuum limit that is consistent 
with the experimental value of the mass ratio. The same holds true for m^/mp. 

In order to compute the masses of strange hadrons as well, Butler et al. assumed first 
order S'f/(3)flavor symmetry breaking. These results, once they are extrapolated to the 
continuum, are in agreement with the experiments as well. 

In all these cases, the error bars are somewhat large and the number of points is too 
small to check whether the mass ratios really depend linearly on the lattice spacing. 
Nonlinear terms could be important in this range of lattice spacings. To some extent, 
we can investigate this question by using another precision-quantity: the potential 
V{r) between static quarks can be calculated very precisely in lattice QCD. A well 
calculable length scale tq may be obtained from the force F{r) = -§;:V{r) through 
the implicit definition F(ro)rQ = 1.65[^. The constant 1.65 is chosen such that one 

^^ There are strong arguments pq| that the quark-mass dependence of hadron masses is more 
comphcated in the quenched approximation. It may be non-analytic at zero quark mass. It is, 
however, not obvious at which values of the quark mass such behavior sets in. Furthermore, such 
non-analyticities would mean that observables that are calculated in the quenched approximation 
deviate very significantly from the true result when the quark masses approach zero. Therefore, one 
should stay away sufficiently far from the chiral point and extrapolate using the mass dependencies 
expected for full QCD. 
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mpa 

Figure 4: Extrapolation of the nucleon to p mass ratio to the continuum hmit. The figure is taken 
from ref. Ml. 



obtains tq ~ 0.5 fm when one uses the force derived from potential models that success- 
fully describe the bb spectra (0.5 fm is about the distance, where the phenomenological 
potential models are most severely restricted by the spectrum), tq has not been de- 
termined at all values of the bare coupling which we need here. We therefore take the 
following approach: at distances r^ > 0.3 the potential is well described by 



V{r) = 7r/(12r) + 



ar 



(2.34) 



with a representing the string tension^. We use a as defined in eq. ( |2.34| ), noting 
that through the very parametrization of the potential one has ^/arQ = 1.18. For the 
precision that is required in the following, the uncertainty due to the parametrization 
eq. (|2.34|) is irrelevant. In the future, when hadronic masses with higher precision are 
available, it will be important to use the direct definition of tq. 



We have combined the data obtained from the simulations of ref. [^, ^Tj as listed in 
p3| with masses of the p-meson from ref. |54, 53, 44, |59||P^ to form the ratio Mp/yT,. 
(By S we denote the estimate of a at a finite value of a in lattice units: S 



ala.: 



^^ The subleading term 7r/(12r) is the universal first correction of an effective bosonic string |5q|. It 
is both in agreement with SU(3) gauge theory simulations and with the force in the continuum limit 
oftheSU(2)-theory|30|. 

^^ In three cases (for L/aMp < 6) it was necessary to correct the mass of the p-meson for a 
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Figure 5: The lattice spa cing dependence of MpjyfT, for three different fermionic actions: Wil- 
son action (full squares ||54 ^ Q, ^), Sheikholeslami Wohlert action (triangles |^, |6^, ^) and 
staggered fermion action |Q (filled squares |6l|, |6^, |6^, |64|). 



The lattice spacing dependence of Mp/yS is shown in fig. |^ (open squares). A hnear 
relation describes the data well in the whole range. The range of lattice spacings that 
is covered here includes the one in the extrapolations of [Q but it extends further 
towards a = 0. The linearity observed in fig. |^ lends support to the extrapolations! 
that we discussed above. 



The extrapolation of Mp/V S is further checked by computations that were done using 
different forms of the fermion action. 1) The p-mass data of various groups |^l], ^, 
63| , |64| using the staggered fermion action [Q (filled squares) have corrections to the 
continuum limit of the opposite sign. For this fermion action, mass ratios should 



approach their continuum limit with corrections proportional to a^|^. This form 



slight finite volume effect | 
account. 



The correction was taken from 1 60 and its uncertainty was taken into 



23 



is fitted to the data. As a result we observe universality in the continuum limit 
despite the quite significant differences at finite values of the lattice spacing! The 
figure underscores the importance of performing such extrapolations in order to obtain 
the physical continuum limit numbers. 2) The recent computations by the UKQCD- 
group[^ and the APE-groupp6|, |67[ , using the improved action are represented by 
the triangles. Although the error bars are somewhat large to draw definite conclusions, 
it appears to be close to the above continuum limit extrapolations indicating small 
discretization errors. Further points for larger values of the lattice spacing in this plot 
could clearly settle the question to which extent the improved action of Sheikholeslami 
and Wohlert reduces the 0(a) scaling violations in nonperturbative quantities. 

Altogether the investigation of the lattice spacing dependence and universality of 
Mp/yTj provides further support for the continuum limit extrapolations of mass ratios 
of Butler et al. As the extrapolated numbers agree with the experimental observa- 
tions to within ~10%, one may conclude that the quenched approximation is at least 
a good model for QCD. In fact, it appears to be the best model of QCD that we have, 
considering that it contains only the fundamental parameters of QCD as unknowns. 
However, we would like to emphasize again that there is an unknown systematic error 
involved when we replace QCD by this model for calculating a new quantity. 

We do not want to review the numerous efforts that have been made to directly 
compare QCD and its quenched approximation (at a given value of the cutoff). Instead, 
we refer the interested reader to review articles on this subject pT|, [5^ . Sea quark 
masses of around the strange quark mass and higher have been used so far. For these 
quark masses, no significant changes in the physical observables have been found apart 
from the aforementioned change in the coupling constant at the cutoff scale. This is 
in agreement with the conclusion drawn in the previous paragraph. 



2.7 The b— Quark on the Lattice 

In the previous sections, we have discussed the general problems in practical lattice 
QCD calculations. Apart from renormalization of operators, the dependence of lattice 
results on the spacing a is most important. As long as we are considering hadrons 
consisting of light quarks such as m, d and s, it is sufficient to have a good resolution of 
the wave-function in order for a-effects to be moderate and to be able to extrapolate 
to the continuum limit. 

When one considers bound states of a heavy quark like the 6-quark with other con- 
stituents, it is necessary to have 

rrih a « 1 (2.35) 

in addition. (We denote a generic heavy quark by h.) In a given background gauge 
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field, the heavy quark only propagates correctly under this condition. 

This condition is hard to fulfill for the following reason. In order to control finite size 
effects in hadron masses and matrix elements, we need to compute with a box size of 
the order of more than one fermi (see section |3.2|) . Even a lattice with L/a = 32 points 
has rrih a > 1 for these values of L. Consequently, a direct simulation of the 6-quark 
requires finer lattices than the ones that have been technologically accessible so far. 
One had to search for alternatives to a simulation with fully propagating 6-quarks. 

Intuitively, one expects that (for low energy matrix elements) the high frequency parts 
of the propagation of a heavy quark can be absorbed into local terms of an effective 



actionpH]. It should therefore be possible to work with a modified action and a cutoff 



a ^ << rrif, as long as one is only interested in long distance observables. 
Through a formal expansion in l/rrih one obtains the non-relativistic QCD action 

Snrqcd = <^^[^o - ^ - ^]0 + Oil/ml), (2.36) 



written here in terms of the covariant derivative D^ = {Dq, D) |2^ and the chromo 
magnetic field strength, B. In this formulation, the heavy quark field is a 2- 
component spinor field and a is a vector composed of the three Pauli matrices. 



Beyond the formal level, the use of the effective action eq. ( 2.36| ) is nontrivial. The the- 



ory defined through Snrqcd and its appropriate regularization is nonrenormalizable. 
A finite cutoff has to be kept. The condition a~^ << ruh requires that computations 
have to be done just in the opposite range from eq. ( |2.35| ). This means the lattice 
spacing is in the middle of the range in fig. ^. For such resolutions, one has to worry 
about the effects of order 0{a) as it is clearly demonstrated in that figure. These 
cannot be removed by extrapolation a ^ since that limit does not exist. Instead, 
0(a)-effects have to be reduced by adding new (higher dimensional) operators to the 
action. Their coefficients need to be determined by matching a number of observables 
in the effective theory to the observables in the full theory. One possibility is to match 
directly to experimental observations. This procedure both reduces the predictabil- 
ity of the theory and introduces statistical errors into the coefficients. The second 
possibility is to estimate the coefficients through mean field theory or by matching 
perturbatively [|69| . 

Concerning the second choice, it is likely that incalculable nonperturbative contribu- 
tions in the coefficients induce finite terms into the final physical matrix elements that 
are of the same order as the physical result that one wants to calculate |]TD[. This 



possibility originates from the powerlaw divergencies of the effective theory |70[]. It is 
argued, however, that such terms should be small numerically |]7l], ^ 



25 



According to our judgment, it will be difficult to quantify the uncertainty due to 
either missing higher order terms in the action or the incalculable nonperturbative 
contributions to the coefficients of the effective action. Currently, non-relativistic 



QCD on the lattice is under intense investigation! 72 



It is important to notice, that these problems of non-relativistic QCD do not exist 
when one truncates the action eq. ( p.36|) after the ffist term|^3|, ^. It has been 
suggested by Eichten|7^ to apply this 0((l/m/i)°) approximation to the b-quark on 
the lattice. This treatment of a heavy quark is called the static approximation, 
since the heavy quark is represented by a static color source |75|, |7^ . 
Sh{x, y) in a given background field is just a Wilson line: 



Its propagator 



Shix,y) = %- (l + 7o)/2 W^{xo,yo;x) for Xq > yo ; (2.37) 



W{xo,yo;x) 



yo-a 



n Uo{z = {Zo,Xi,X2,X3)) . 



zo=xo 



In the static approximation, the only powerlaw divergent - and therefore incalculable 
in perturbation theory - renormalization is the renormalization of the quark mass. 
The latter is of limited physical interest. Most importantly, the renormalization of 
the axial vector current can be calculated in perturbation theory []73|, |7^ allowing for 
a computation of the decay constant in this approximation. 



One can therefore compute properties of B-mesons through the following strategy. 
As a ffist step, one determines their limiting behavior for mh -^ oo using the static 
approximation. Then one investigates their dependence on the mass of the heavy quark 
for masses as large as possible and finally one interpolates between the results at finite 
values of the meson mass and in the static approximation. The interpolation is done in 
the form suggested by the 1/m/j-expansion, but the coefficients of the different powers 
of l/rrih are not computed explicitly since they suffer from powerlaw divergencies [|7D|. 
Rather they are taken from the phenomenological matching between finite mass and 
infinite mass results. 



26 



3. The Leptonic Decay Constants /o, / 



Ds 



This section is a review of the computations of heavy hght decay constants for heavy 
quark masses around the mass of the c-quark. As a general orientation, we hst in 
table m the results that have been quoted in the literature. Note that the different 
computations have used different renormalization constants and different ways to de- 
termine the lattice spacing. At finite values of the lattice spacing and in the quenched 
approximation this can introduce quite significant variations. This explains the spread 
of the values given in table |I[ 



Ref. 



/D[MeV] 



JdsI f. 



D 



a-'[GeV] 



m 



m 
m 



m 



m 



m 



n 



44 



174 ± 26 ± 46 

197 ±14 
181 ±27 
190 ±45 

198 ±17 
210 ±15 

208 ± 9 ± 35 ± 12 
185: 



^42 
-7 



170 ± 30 



(/d, = 234 ± 46 ± 55 MeV) 



1.17 ±0.22 
(/^^ = 209 ± 18 MeV) 

1.08 ±0.02 
1.11 ±0.05 

1.18 ±0.02 
1.09 ±0.02 ±0.05 



2.2 

1.9 

2.6 

1.9 

1.9 
1.9-3.4 

2.9 

1.9-2.6 

1.1 - 2.8 



Ui 














m 



283 ± 28 
130 - 300 



(/^^ = 170 - 395 MeV) 



1.5 
1.5 



Table 1 ; A compilation of published values of the decay constants in the D-system. Note that the 
treatment of systematic errors is very different in the different works( In particular, the large range 
quoted in [Q is due to a conservative estimate of the systematic errors, while [Q did not estimate 
systematic errors.) All computations use the Wilson action {csw = 0) except for ref. |46 , which 
uses the treelevel 0{a) improved action(csiy — 1)||2J. The forth column gives the value of the cutoff 
estimated from the string tension as a^^[GeV] = 0.420 /vS . 



In order to combine the information from the different publications, it is necessary to 
separately investigate each of the possible sources of systematic errors. A good statis- 
tical accuracy is needed before systematic errors can be investigated in a meaningful 
way. The pioneering computations |7^, |7D|, ^, |T^ did not have a sufficient accuracy 
for that purpose. Therefore, we will only use the results of refs. ||8^, ^, ^, |7^, ^ in 
the following. It also would have been of great use to include the computation that 
was performed at the smallest value of the lattice spacing[^. Unfortunately, applying 



what we learned in section |2.5| about the effects of excited states to that computation 
shows that systematic errors due to excited states are of the order of the statistical 
errors. Consequently it is not included in the following. The same holds true for part 
of the results of ref. |8^; the conclusions of that paper have to be revised. The data 
that we do include in the following discussion fS^, |46|, 



|, IJ, |73] either was obtained by 
appropriate smearing techniques or at sufficiently large separations Xq. Uncertainties 
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due to excited state contributions should be well below the statistical errors. 



3.1 Extrapolation to Light Quark Masses 

In the work mentioned above, all extrapolations are performed as follows. Motivated 
by chiral perturbation theory P3 one fits 



MUgl Kf, Kj,) = A{gl){-^ + ^--J-) (3-1) 

ZKf ZKf> Kcrit 

to determine A{g^) and Kcrit- Then the points on the renormalization group trajectory 
corresponding to the s- quark and the u/d-quark may be fixed by requiring the pion 
and the Kaon to acquire their physical mass in units of the string tension or another 
dimensionful observable such as /^r 

A{9l){^--^)/^ = ml/a (3.2) 

with m^ = 137 MeV, tuk = 490 MeV and a ~ 420 MeV. 

Other quantities, such as masses and decay constants of HL-mesons, have a much 
weaker dependence on the bare quark mass. They are extrapolated linearly in 1/Kf to 
the corresponding points Kf = Ku or Kf = Kg. The linear behavior is confirmed in all 
simulations (within the errors). So, for the quantities discussed here, the extrapolations 
appear to be quite reliable and are not discussed any further^. 

We proceed now to investigate finite size effects, the lattice spacing dependence and 
renormalization in the quenched approximation to arrive at an estimate of the decay 
constant in that model. In a final subsection, we will try to see to what extent light 
sea quarks alter that result. 



3.2 Finite Size Effects 

At asymptotically large lengths L of the box, the dominant finite size effects origi- 
nate from a self interaction of the meson through an exchange of a particle with the 
quantum numbers of the vacuum. The particle can propagate once around the world, 
introducing an L-dependence. In the quenched approximation the relevant lightest 



^"^ We note that the choice of scale (namely a in the above equations) influences the precise values 
of Ku and Kg. Numerically, it is not important in the following, however. 
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particle is the 0"*"+ glueball. The effect is suppressed exponentially in L with a decay 
rate of the inverse glueball mass. For lattice sizes that are being used in the simula- 
tions of HL-mesons the product LmQ++ is larger than ten[^, |8^ 
finite size effects is therefore irrelevant in the present context. 



This mechanism for 



The dominant finite size effects within an intermediate regime should originate from 
a distortion of the wave function in the finite volume. In a nonrelativistic potential 
model, the rms-radius of an HL-meson is significantly smaller than the rms-radius of 
a LL-meson since the reduced mass is smaller by almost a factor of two. Taking a look 
at LL-meson decay constants that are published in the literature, a size-dependence 
is not visible for LyS ~ 3 and larger. As a result, we clearly do not expect relevant 
finite size effects for HL-mesons once L is in that range. 
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Figure 6: The dependence of the HL pseudo scalar decay constant F on the box size LQ. The 
hght-quark mass is fixed at about twice the strange quark mass and the meson mass varies from 
about 0.8 GeV (circles) to 1.5 GeV (filled boxes). 

Indeed, as shown in fig. P, numerical investigations with moderate light-quark masses 
and varying mass of the heavy quark, confirm the expectation E^ . 
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One may speculate further that the finite size effects seen within the regime of LyU ~ 
1.5 — 2 originate from a distortion of the hght-quark wave function in the finite 
volume. For an exponentially decaying wave function, the finite size effects on the 
wave function at the origin are again exponential in L. The characteristic scale is the 
coefficient in the falloff of the wave function. A rough estimate of 1.5 in units of the 
string tension is derived from the Bethe-Salpeter wave function of HL and static-light 
mesons in Landau gauge |^. The data in fig. ^ suggests an even larger coefficient of 



the falloff if one fits to an exponential law in the whole range. Consequently, finite 
size effects are expected to be smaller than the statistical errors when La/E ~ 3. They 
are neglected in the following. This statement is checked by the data on the level of 



5%[^. Note that finite size effects must be reevaluated once future investigations 



reach higher precisions. 

3.3 Lattice Spacing Dependence and Renormalization 

The most important systematic uncertainty of the decay constant originates from the 
dependence on the lattice spacing. In principle, the cleanest way of investigating this 
dependence and altogether removing the uncertainty through an extrapolation to the 
continuum is the following. 

One fixes the hopping parameter of the light quark Ku as described above and the 
hopping parameter of the charm quark Kc is determined through 

Mp{Kc, Ku, gl)/Fp{Ku, Ku, gl) = nio/U . (3.4) 

At these values of the hopping parameters one considers the ratio 

Qigl) = Fp{Kc, Ku, gl)/Fp{Ku, K^, gl) . (3.5) 

Up to lattice artifacts, this is the ratio of the D-meson decay constant to the pion 
decay constant f.,,: 

Q{gl)*U = fD + 0{a). (3.6) 

In the literature this is called "setting the scale through /^r". The advantage over 
every other method is that the renormalization of the axial vector current (see section 
T^ cancels, apart from its dependence on the mass in lattice units, which represents 



a lattice artifact. This means that the axial vector current is renormalized nonpertur- 
batively. 

With increased statistical precision, this method should give the most reliable deter- 
mination of fo- At present, the error bars of Fp[Ku, Ku, g^) are relatively large and 
dominate the errors of the ratio Q. Therefore, the 0(a) terms cannot be detected well 



enough to allow for a continuum extrapolation using eq. (|3.6|) 
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Consequently, it is advantageous to follow a detour when one wants to exploit the 
results presently available. We split Q into 

Q{9l) = qi{9l)/q2{gl), (3.7) 

q,{gl) = Fp{Kc,K^,gl)/Vf., q^igl) = Fp{K^,K^,gl)/Vf. . (3.8) 

Then qii^g^) and q2{go) can be extrapolated separately to the continuum limit and fa 
is calculated by their ratio in the continuum limit. 

Since the renormalization of the axial vector current is known to one-loop order only, 
there are in addition to the terms of order 0{a), errors of order O(a^) in gi and q2- 
From the work of Lepage and Mackenzie [p5|] , it appears justified to assume that these 
missing terms have coefficients of order one. A coefficient of four would mean 4 a^ ~ 
2% to 8% in the range of cutoffs that are considered below. As it will become evident 
below, an uncertainty of that order does not need to be too disquieting. In addition, 
the same systematic uncertainty is present in gi and q2 and the effect of the uncertainty 
on the extrapolation will cancel partly when one takes the ratio. 

The two advantages of this procedure are that i) the lattice spacing dependence of gi 
can be studied with small statistical errors and ii) additional data at other values of g^ 
are available for the pion decay constant. Due to this additional data, most notably 
the one of ref. |8^, g2 can be extrapolated quite well. To do this, we have taken 



the data for /^ from the literature [^, ^, |^, |^, ^] and applied the perturbative 



renormalization of the axial vector current eq. ( 2.24 ) 



We show the extrapolation of g2 to the continuum in fig. |^. The lattice spacing 
dependence is very weak and we obtain in the continuum limit of the quenched ap- 
proximation: fn/y/o' = 0.270(12). The error does not contain the (maybe around 5%) 
uncertainty due to the missing terms of order O(a^) in the perturbative renormaliza- 
tion. 

Fig. 1^ also contains results obtained with the improved action. Within their uncer- 
tainty they coincide with the Wilson-data. This is in agreement with the assumption 
that both lattice spacing effects and higher order perturbative terms in the renormal- 
ization are small. 

The extrapolation of gi is a more difficult task. Fig. ^ summarizes the available data 
of sufficient precision. The filled squares are results with the Wilson action, the nor- 



malization eq. ( p.21|) and 1-loop perturbative Za- The open squares are the same 
data with normalization eq. ( |2.22| ). It is more than apparent that the lattice spacing 
dependence is stronger for this normalization. This matches completely with the case 
of the local vector current renormalization discussed in the appendix. From a theoret- 
ical point of view, one can only speculate about this surprising result. It is helpful to 
recall that lattice artifacts originate both from the operators in a correlation function 
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Figure 7: The lattice spacing dependence of the pion decay constant. The data for the pion decay 
constant was obtained from several simulations. For the Wilson fermion action (squares): enumerating 
the points from left to right, ref. |82 contributes to points 1,5,8; ref. |Q to 4,6,8; ref. Q] to 2,5,7 
and ref. [p3[ to point 5. Point number 3 as well as the points with the SW action (triangles) are taken 
from ref. |l59|. 



and from the correction terms in the action. While eq. ( p.22| ) is designed to improve 
the operator at treelevel, the action remains unchanged. A possible interpretation 
of fig. 1^ is that without this "improvement" the lattice artifacts originating from the 
operator and the one from the action partly compensate each other, such that the 
full squares show little a-dependence. On the other hand, the open squares show the 
lattice artifacts originating from the action. In this interpretation the latter show the 
"natural" size of lattice artifacts which appears to be somewhat reduced when one 
systematically improves both the action and the operators (triangles). 



Without performing any fit, by mere inspection of the data points, all three data sets 
are consistent with extrapolating to the same value in the continuum limit S = 0. The 
question is, whether any one of the data sets can be extrapolated to S = reliably. The 
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Figure 8: The lattice spacing dependence of the D-meson decay constant. The data with the Wilson 
action are from ref. |Q (first point) ref. [B^] (third point) and ref. [Q. The triangles represent the 
results with the improved action. They are from ref. pTj (first point) and |46[ (second point). 



large difference of the different data sets simply reflects that the mass of the charm 
quark in lattice units is quite large at least for \/T, > 0.3. Therefore, we omit the 
point with the largest value of the lattice spacing and extrapolate the others linearly 
in vS- The points at S = in fig. ^ show the results of the extrapolations. They 
are consistent and furthermore, results within those error bars would be obtained if 
one included the points at the largest value of the lattice spacing. We take this as 
evidence that the error bars include the systematic error from the extrapolation. In 
particular, the full squares extrapolate to /d/v^ = 0.338(39). The uncertainties due 
to the perturbative renormalization are as above. 
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3.4 Estimate in the Quenched Approximation 

We combine /d/v^ = 0.338(39) and /^/v^ = 0.270(12) to attain 

Jd/U = 1.25(15) . (3.9) 

Here, we have neglected the O(a^) renornialization error. This uncertainty, which may 
be sizeable in the individual ratios qi, should cancel out to a large degree when one 
forms fn/ fiT- In particular, we may also regard the fitted function in fig. |^ as a legiti- 
mate interpolation of the data. Since it does not have any significant a-dependence, 
the ratio Fo/vS contains to a reasonable approximation the full a-dependence of 
the ratio Q. As explained above, Q has only 0(a) effects and no 0{g'^) corrections. 
Consequently, the linear extrapolation of gi and g2 separately seems justified at the 
present level of precision. 

The above value corresponds to fn = 164(20) MeV. 

Regrettably, the computations with the improved action do not yet cover a sufficient 
range of lattice spacings to perform such an analysis. As seen in figs. ^ |^, the data 
are in agreement with the ones obtained with the Wilson action so far. 



3.5 Full QCD 

Simulations of full QCD have so far mainly been carried out using the staggered 



fermion action||9^. In this case, one has fewer degrees of freedom, which renders the 
simulations more tractable. Unfortunately, no results for HL-mesons exist with that 
discretization of QCD. It appears to be particularly difficult to find operators that 
have sufficient projection onto the ground state in the HL channel when using this 
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In two publications, listed as the last two lines in table |l], decay constants of HL- 
mesons are computed with the Wilson action. Ref. |]S3| presents a very exploratory 
simulation with an approximate algorithm to account for the fermion determinant. 
On the other hand, ref. |]7^ simulates with the Hybrid Monte Carlo algorithm with 
two different sea quark masses, one of which is some 25% below the strange quark 
mass and the second one is ~15% above the strange quark mass. The lattice spacing 
determined from the mass of the p-meson is around a = 1/(1.5 GeV). This corresponds 
to v^ ~ 0.38, the point with the largest lattice spacing in fig. ||. Consequently, the 
systematic errors due to a finite lattice spacing are estimated to be rather large ||77||, 
which results in the wide range for fjy, fn^ that is quoted in table |I|. 

We want to go somewhat beyond the interpretation of the data presented in ref. 
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Figure 9: 



The ratio Fp{gQ,Ki,Kh)/Fp{gQ,Ku,Ku) as a function of 



Mv{gQ, Ku, Ku)/Mp{gQ, Ki, Kh). Open symbols are for nj = 2 |r^ with Ki somewhat below Kg 
(squares) and Ki somewhat above Kg (triangles) . Filled circles are for nj — G and Ki = Kg [Q . 



For the hadron mass spectrum, dynamical fermion simulations have yielded results 
that are completely compatible with the quenched spectrum at the same value of the 
lattice spacing. It is in fact quite plausible that one should compare at a fixed value 
of the cutoff, since if one has a small effect of quark loops on a certain observable it 
is natural that also the lattice artifacts do not feel much of the presence of dynamical 
fermions. 



In order to test to what extent this statement is correct for HL decay constants, we 



compare the numbers given in the tables in ref. [^ with the quenched data of ref. [Q 
at a/S ~ 0.38. The systematic uncertainty on Za is eliminated by looking at the ratio 
of HL-decay constants to the pion decay constant]^ A further uncertainty originates 

^^ We used the normalization eq. ( ^.2l| ) . Inserting eq. ( ^.22[ ) would result overall in quite different 
numbers, but the agreement of the quenched QCD with the full QCD results would be just as good. 
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from the masses of the quarks. In fig. |^, we show the quenched result for the strange 
quark mass and the results with Uf = 2 for the two sea quark masses that straddle 
that mass. The heavy quark mass varies along the x-axis where the mass of the meson 
is given in units of rrip. 

Remembering that the quenched mesons have a light-quark mass in between the two 
(independent) simulations with dynamical fermions, there is complete agreement be- 
tween the two sets of data. We may take this as positive evidence - albeit with a 
precision of maybe 10% - that there is no effect of dynamical fermions on the decay 
constants. Of course, this is a check at one relatively large value of the lattice spacing 
and one cannot exclude a somewhat different result at a smaller value of a. We should 
remember further that the light quarks have masses around the strange quark mass in 
this test. 

Further tests with smaller a and smaller quark masses will become possible soon. 
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4. The Leptonic Decay Constants /g, / 



Bs 



Computations of properties of B-mesons are considerably more difficult than the ones 
discussed in the previous section because of the large mass of the b-quark (cf. section 



2.7). At present, appropriate values of the lattice spacing, i.e. such that rriha « 1, 



have not been reached. 

Therefore, the estimation of fs starts from the limiting value obtained in the ap- 
proximation of neglecting terms of order 0{m^^). Next, one investigates the mass 
dependence of HL decay constants for varying masses of the heavy quark, where those 
masses are around the mass of the charm quark. Finally, one matches the two results 
through some interpolation and arrives at an estimate of Jb- We describe these steps 
in the following three sections. 



4.1 Static Approximation 



Using a static propagator eq.( |2.37| ) for a B-meson, the (bare) axial vector correlation 
function does not depend on the mass of the heavy quark. So, the combination Fp\/Mp 
does not depend on the mass of the quark in the bare regularized theory. Since we are 
working with an effective action for the heavy quark, we have to restore the relation 
to QCD by matching the effective theory to QCD [|73|, IT^I- It is most transparent to 



require directly that the correlation function eq. ( p.l3| ) is the same in the effective 



theory and in QCD up to terms of order O^mJ^^) |7^. To 1-loop order perturbation 
theory this is achieved after accounting for a (linearly divergent) subtraction of the 
binding energy and a renormalization of the axial vector current by a factor ||73|, 173, 0^ 



■^statl*^ ^h) = 1 + a [-log(a nih) - 2.372] ,for csw = , (4.1) 

vr 

^stat(« mh) = l + a [- log(a rrih) - 1.808] , for csw = 1 • (4.2) 

n 

In the static approximation one thus computes 

F^^""^ = ^stat(« ^b) < 0|g,7o75gz(0)|P > |stat , (4-3) 

where, because of the physical interest, we have renormalized at mi,, p contains the 
full information on the decay constant that can be obtained in the static approximation 



in the formPS 



F^t^t = lim F{Mp) (4.4) 

A/p— »oo 

F{Mp) = P(Mp)^(^^^^ °; Po = II -luf, (4.5) 

37 



where nj is the number of dynamical quark flavors (Here the leading logarithms have 
been summed using the renormalization group equation and the difference between 
the mass of the meson and the mass of the quark has been neglected) . 

As a first overview of the effort that went into the computation of F , we list in 
table the different published computations of f^^^^ = a^^'^ F / y/rriB. Here, one 
has to note, that part of the scatter in the table originates from different ways of 
setting the scale and other systematic errors, which we attempt to analyze below. 



Ref. 


/gat[MeV 


fst&t / -fstat 


a-^[GeV 


csw 


1941 


<~260 




1.9 





18 


310 ± 25 ± 50 




1.9 





^ 


366 ± 22 ± 55 


1.10 


1.9 





43 


230 ± 22 ± 26 


1.16 


1.4-2.7 





66 


350 ± 40 ± 30 


1.14 ±0.04 


1.9 






95 



82 



66 



H 



319 ± Xp^g X iij^Q^y, 
235 ± 20 ± 21 
370 ± 40 
253: 



3/2 



1.11 ±0.04 



-16 
-15 



-105 
-14 



1.14: 



-0.04 
-0.03 



Table 2: A compilation of published values of the decay constants in the B-system as they were 
obtained in the static approximation. The fourth column gives the value of the cutoff estimated from 
the string tension as a~'^[GeV] = 0.420 /V^ ■ 



Ground state domination 

The question of how to extract the ground state meson decay constant has been dis- 
cussed quite extensively over the last years. The problem is described in section p75 . 
Here, we summarize what is known in particular for mesons in the static approxima- 
tion. 



Following the original suggestion of Eichten ||75[ to compute F , it was noted by 
Boucaud et al. |Q, that when one employs the static approximation, the statistical 
fluctuations increase rapidly with the euclidean time separation xq- Therefore, with 
local operators one never observes a signiflcant plateau in the effective mass. For that 
reason, they could only give an upper bound for the decay constant. 



Subsequently, nontrivial wave functions were used||45 
be estimated in 



1^ 



T|, ID and F^^^^ could flrst 
1^. Quite some effort went into the construction of good wave 
functions for static-light mesons. Nevertheless, in all computations that have been 
done so far, one sees only relatively short plateaus in the effective masses. Ref. [Q 
stimulated the discussion on this point. At present the best solution to the problem 
is given in ref. [95]. There, the smearing wave function is determined through the 
diagonalization of a matrix correlation function. One thus obtains an optimal wave 
function out of a space of linear combinations of a number of basis functions. As a 
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result, out of all computations, the plateaus in the effective masses are longest and 
most convincing in ref . |^ . Unfortunately, only a result at one value of a is presently 
published with that method. 

A point of discussion is also how to fit the two correlation functions C^'^ and C^'^°'^ in 
order to extract the matrix element of the local axial vector current. The analysis of 
refs. [^, ^ is based on a combined fit to C^'^ and C^'^°'^. The appropriate range of the 
fits is determined from their x^. Since a large correlation matrix has to be determined 
from a quite limited statistical ensemble, this criterion may be misleading! 



97. 



The other groups determine the (linearly divergent) binding energy Mp from the cor- 
relators C^'^, where / labels the smearing wave function. These correlation functions 
have a positive spectral representation and the corresponding effective mass is easy 
to interpret. Also, because it is the correlation of two smeared fields, the unwanted 
excited states are more strongly suppressed than in C^'^°'^. The analysis then proceeds 
in the following way: either one uses the binding energy as extracted from C^'^ , per- 
forms a fit to C^'^°'^ constraint by this energy and then obtains the decay constant 
from the amplitudes in the two fits|T9| or one determines the height of the plateau 
of the ratio c^'^°'^ /C^'^ and combines that with the amplitude of the fit to C^'^ . For 
further variations on the theme, see ref. |]43| , Q. Concerning this analysis method, 
one may argue [B2] that the determination of the binding energy from C^'^ is somewhat 
uncertain because there are larger statistical errors in C^'^ than in C^'^""^. 



Altogether, some amount of ambiguity remains and one may suspect that some of the 
published results have an underlying systematic error which might be even as large as 
the statistical error. 



We note, however, that there are nontrivial checks on the results 
(3 = 6.0, csw = 0, there are five computations 111, [ll, |6^, [ 



In particular 
H. Ref. |6 



at 
finds that 



the fits in |T^ started somewhat early, introducing a small systematic error. Ref. [Q 
extends the calculation of ref. |T^. Therefore, we compare []66| , |43| , |8^ . The interesting 
point is that three different types of wave functions were used. Refs. [^ ^ used a 
wave function which is a three-dimensional cube for the light quark, implemented in 
Coulomb gauge. The sizes of the cubes were 5-7 lattice spacings|]BB| and 9 lattice 
spacings|8^. In obvious notation we label them by C5 — C9 . Ref. |^ used two 
different, gauge covariant, wave functions. One has an approximately gaussian shape 
('G') and the other one an approximately exponential shape ('-£")• ^^ can compare 
the results from the three different wave functions directly at one value of the light 
quark hopping parameter K = 0.1540. The results quoted are < 0|g/j7o75g/(0)|P >g= 
0.36(3), < 0|g;,7o75g/(0)|P >e= 0.37(2) H and < 0\qhlol5qimP >C5= 0.41(5), 
< 0|g/i7o75g/(0)|P >C7= 0.38(3) ||HB[. So the decay constants are in good agreement 
for these four different wave functions. Ref. ^2[ obtains a somewhat low value < 



0|g/i7o75gz(0)|P >C9= 0.31(4) due to the analysis method discussed above. Moreover, 
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G.loci 



Xo) 



it is interesting to note that the systematic effect that arises from a fit to C 
at too small a value of Xq is to decrease the value of the decay constant, while for 
C^''°^(xo) this results in an overestimate of the decay constant. So the agreement 
that is found at large values of xq is nontrivial and suggests that xq is large enough in 



the analysis of ref. ||^. Also in the computations with the action of Sheikholeslami 



and Wohlert [^, |66[| , two different wave functions are seen to give consistent results at 
/? = 6.0. These consistency checks indicate that the problem of contaminations from 
excited states is not too severe. 
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Figure 10: J^^tat -^^ units of the string tension as a function of the lattice spacing. Data of 
refs. p3| , p3, pS, p2l M are shown and have been combined, when they were available at overlapping 
values of the lattice spacing. Open boxes correspond to the Wilson action and filled triangles to the 
improved action. 



Finite size effects of F were investigated in refs. |^, |^ . We can already infer 
from fig. H, that they should be negligible once L-y/S = 3 (the figure shows that the 
finite size effects decrease as the heavy quark mass increases). This is in agreement 



with the study of refs. [^, p5[] . 
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Lattice spacing effects are a significant remaining source of a systematic error. In 
order to investigate it, we consider _pstat ^^ units of tlie string tension. The data 
of refs. p3| , p5| , [66| , ^, ^ is shown in fig. [1^. Here, the renormahzation constant 



was obtained by using the perturbative expansion in terms of a. Overall, the lattice 
spacing dependence appears not to be very strong. This statement refers precisely to 
the ratio j^stat^^3/2. j^^ most of the early computations, the mass of the p-meson was 
used to set the scale. Due to the strong dependence of Mp/\/E (with csw = 0) on 
the lattice spacing, there is a strong dependence on a when the scale is set through 
TUp. This explains why the earlier computations, working at vS > 0.2 and setting the 
scale through rrip, obtained quite large values of F^^at ^^^^ table ^. 

The data with csw = is extrapolated to zero lattice spacing as shown by the dashed 
line: y'Stat^^3/4 _ i,35(^i7^_ One can see at a glance that the fit does not have a good 
X^- However, it is obvious from the figure that this is not due to deviations from the 
assumed linear dependence on the lattice spacing. Rather, it is due to the scatter of 
the data itself. The most likely explanation for the scatter is incomplete ground state 
domination in some of the data points. A systematic error of the order of half of the 
statistical errors could account for the scatter in fig. ITU. 



Since a better procedure is lacking at the moment, we account for this scatter by 
taking in addition to the statistical error of the extrapolated value a systematic error 
of the same size. In the future, this systematic error should be removed by a further 
improvement of the trial wave functions. 

We further estimate the error of the missing higher order terms in the perturbative 
expansion of Zgtat in the following way. We allow for a 2-loop term 5 a in Zstat 
(a coefficient as large as five is certainly a conservative estimate). With this term 
the whole extrapolation is repeated. The change in the extrapolated value gives the 
systematic error ±0.08. 

The data with the improved action {csw = l^ triangles) are fully consistent with the 
Csw = points. We must, however, remember that in the comparison between the 
two sets of data an 0(Q;^)-uncertainty exists and the a-dependence could be quite 
different. We should only compare the extrapolated values. With the two points at 
hand, such an extrapolation is not possible]^. 

Altogether we obtain 

jstat/^3/4 ^ 1.85(17) (17) (8) . (4.6) 

As in the case of fo, it is better to change from a prediction in units of the string 
tension to one in units of /vr, since a phenomenological value for the string tension is 



^^ Clearly, there is no evidence that the a-dependence in this observable is reduced compared to 
csw = 0. 
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based on assumptions. Furthermore, one may hope that effects of dynamical fermions 
are reduced when one normahzes to a decay constant. So we combine y'Stat^^3/4 ^[f^]^ 
/^/y^ = 0.270(12) to 

xstat 

^^7^ X (0.132 GeV)3/2 = 0.63(6) (6) (3)GeV^/2 . (4.7) 

This corresponds to fs = 276(26)(24)(12) MeV + 0(171^/1713). The errors quoted 
are (in order of appearance) the statistical errors including the extrapolation, the 
systematic error due to possible contaminations by excited states, and a 5a^ term in 
the renormalization. (For the latter, it is not justified to assume a cancellation between 
the higher orders in Z^tat and Za because these renormalizations are of quite different 
origin.) 

We may also change to physical units through the mass of the p-meson. Here, the 
extrapolation of fig. || resulted in nip/ y/a = 1.81(4). Therefore, 

fstat 

^^-^ X (0.77 GeV)^/^ = 0.51(5)(5)(2)GeV^/2 (4.8) 

rup' 

is a further estimate corresponding to /^ = 224(22)(21)(10) MeV + 0{mj^/mB)- 



The difference between these two estimates is most likely an effect of the quenched 
approximation||89|. Consequently, the difference between eq. (|4.7|) and eq. ( [4. 81) should 



once again be regarded as a systematic error. 



4.2 Mass Dependence 



In order to get an estimate of the terms of order Olm^/mB), one may examine the 
mass dependence of the decay constant in the mass region around the mass of the 



D-mesonp3|, |8T| , |82| , |46| , |4^ , p% . We do not consider the results of [^, |81| further at 
this point, because they were (partly) influenced by data with contaminations from 
excited states. In addition, the subsequent simulations are more precise. 



Ref. 1^ uses the normalization eqs. ( |2.22| , ^l25|) for the axial vector current. At a fixed 



value of the lattice spacing, they observe that F{Mp) for mp = Mp/a = IGeV - 3GeV 
joins smoothly with p obtained at that value of the cutoff. This is interpreted 
as showing that this normalization is the correct one in the sense of having small a- 
effects. This behavior at a fixed value of a is, however, automatic and is not necessarily 
connected to a reduction of a-effects for intermediate values of the mass. Rather, with 
this normalization, the expression for F(Mp) becomes identical to F if one takes 
the mass of the quark in lattice units to infinity (apart from the renormalization 
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constants that differ slightly). Therefore, F{Mp) automatically joins smoothly with 
the static point. Fig. |^ shows that strong a-effects are present when this normalization 
is used. These effects are also visible in fig. 5.18 of ref. |p2[. We conclude that the 



belief that an application of eqs. ( p.22| , pl25D essentially eliminates the dependence on 
a in the difficult region where ruh becomes close to a~^ is not justified. 



Instead, one needs to extrapolate [^ to a = at a fixed value of mp as it is shown 
in fig. |] for the case of the D-meson. This both requires systematic calculations at 
different values of the bare coupling and increases the statistical error of the final 
result. However, the latter is necessary, since it is only in this way that the errors due 
to the discretization are covered. Comparing different normalization prescriptions or 
results from different actions is indicative at most. It can not replace the systematic 
study of the a-dependence with one action and one normalization condition. 



We have supplemented the extrapolations of [Q with the raw data of |]82| , [67[| . The 



extrapolations are illustrated in fig. ^ for the largest and the smallest value of the 
mass of the heavy quark that was used. All that was said in section ^^ remains true 
for these other values of the mass. We only include the points with a/S ~ 0.4 at the 
smallest values of the mass. This corresponds to a cut on the mass of the heavy quark 
in lattice units of about niha < 3/4. 
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Figure 11: The lattice spacing dependence of HL-decay constants. Data points as in fig. H. 
Changing again to a prediction in terms of /^, we finally arrive at the mass-dependence 
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of fijnp) depicted in fig. ^ 
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Figure 12: / as a function of 1/mp. The two points at 1/mp = represent eq. (4.7) and eq. (4 



All data has been extrapolated to the continuum allowing for a linear dependence on the lattice 
spacing. 

In principle, we would now like to fit the results at finite values of nip together with 
the result in the static approximation to a power series in 1/mp. Unfortunately, we 



have significant systematic errors (included in fig. |T2[), in particular for the result 
in the static approximation. Thus, it is not appropriate to perform such a fit. We 
just connect the upper and the lower ends of the error bars with functions f{mp) = 



Co + ci rrip^ + C2 nip^ instead. 



At mp = mp we read off from the error band 

fp = /b, = 180(46) MeV 



(4.9) 
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One may compare this result to[^ fs/f-K x 132MeV = 186l2iMeV, a value obtained 



with the improved action at one value of the lattice spacing a^^ = 2.7GeV. 

The value quoted here for Jb is obtained in the quenched approximation. The discus- 
sion in section Ol suggests that a similar value would result in full QCD. 



Finally, the columns 3 of table |1] and table ^ may be summarized to 

fBjfB, = 1-10-1.18. (4.10) 

4.4 Comparison to QCD Sum Rule Estimates 

There are a number of QCD sum rule (SR) computations of HL decay constants. 
These results are summarized in table ^ and are compared to the lattice gauge theory 
results in the quenched approximation. 

Ref _ __ Method fn/U fs/U fS'^'/U 

m, m, llOOli Laplace SR 



[ |IOT| , |9| Hilbert Moments SR 

ig, |21|] Laplace SR in HQET 

eqs. ( |3.9| , [4.7| , ^^ LGT quenched approximation 



1.1 - 


- 1.5 


0.9- 


- 1.5 






1.5- 


- 1.9 


1.1- 


- 1.6 










1.1- 


- 1.9 


1.5- 


2.3 


1.1- 


-1.4 


1.0- 


- 1.7 


1.7- 


-2.5 



Table 3: Estimates for HL decay constants. 

We do not want to discuss details of the sum rule approach here. The interested reader 
may, for instance, consult the recent review ref. [|102|] . We point out, however, that the 
error estimates in table ^ are obtained from a variation of the input parameters. They 
do not include uncertainties due to the truncation of the operator product expansion 
or the question of the applicability of perturbation theory in the regime where it is 
used. The latter question is a particularly relevant one, since in the limit of a large 
heavy quark mass, the only relevant dynamical scale is the r.m.s.-radius of the light 
quark wave function. At such a low energy scale, it is very questionable whether 
perturbation theory can be applied. 

An essential difference between the QCDSR estimates and the lattice gauge theory 
approach is that the latter can be systematically improved. We elaborate further on 
this point in the following. 



4.5 Potential Improvements 

The present best estimates for j'b still contain a significant uncertainty. As we have 
explained in the introduction, the value of Jb is very important in the phenomenology 
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of the CKM-matrix. A higher precision is clearly desirable. 

To this end, it is important to note that the overall computational effort that went 



into eq. ( [4 .91 ) is quite small, when one takes a few months on todays massively parallel 
computers as a scale. Nevertheless, the investigations that have been done, have lead to 
eq. ( [4.9|) and, even more important, a semiquantitative understanding of the different 
sources of systematic errors (it goes without saying that the mean values have changed 
due to this). We therefore have a reasonable understanding of how to reduce the errors 
further. 

In order to estimate what can be done realistically, let us assume that a parallel 
computer performing at 20 Gflops is accessible to a group of physicists for a good 
fraction of the time. 

On such a computer with the existing algorithms, the computation of all correlation 
functions on one configuration of 32^ x 96 points at vS = 0.1 takes roughly 6 hours. 
Therefore, one can obtain data points for figs. ^, |TT] which have the average precision 
seen in the plots and which are at a 30-40% smaller value of the lattice spacing. This 



enables one to add points at masses as high as about 4 GeV in fig. ^ with the present 
precision. At the same time, one can reduce the errors of the points in the range up to 
4 GeV because a new point, closer to the continuum, is added in each extrapolation!^. 
In addition, the values in the static approximation can be significantly improved by 
reducing the uncertainties due to contributions by excited states. The latter point is 
already being addressed at Fermilab|103|. 



Combining this information, an overall reduction of the error of Jb (within the quenched 
approximation) by a factor of 2 seems plausible within a period of a year or two. 



^'^ It would also be useful to add points at Ve ^ 0.3 to study the a-dependence further. 
Note that we assume that the computations are performed at yS L ~ 3. With increased precision, 
the size of finite size effects must be checked, of course. This does not pose a problem since it can be 
done at a relatively large value of the lattice spacing. 
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5. B— Parameter 

In the Standard Model, the mixing of Bq and Bq is mediated by box diagrams with 
two M/^-boson exchanges. This corresponds to a short distance operator O^^^^, which 
in the continuum and at distance scales that are large compared to the interaction 
range of the weak interaction has the form 

0^'='{x) = {d{xh,{l - ^,)b{x)){d{xh,{l - ^,)b{x)) . (5.1) 

Its matrix element is conventionally denoted by 

B^^ = a,i^^)-'/''Bsi^^) (5.2) 

with 

^5B(/x)/|M|=<i?o|0^''=2(0)|5o> , (5.3) 

and /i denoting the renormalization scale. In the lattice regularization (Wilson action), 
the operator Of^^^'^{0) is given to lowest order in a by the naive operator eq. (|5.1| ). To 
first order in a, it receives admixtures from operators of the same flavor but different 



chiral structure ||104|| . This mixing is possible because of the breaking of chiral symme- 
try in the Wilson formulation and has been taken into account to first order in a so 
far. 

The B-parameter eq. ( ^Tsf ) at scale /x = a~^ can be computed directly from the ratio 
of three-point and two-point correlation functions 



BBifi = a-'), (5.4) 



E,Ey<M{x)Otr\o)M^{y)> 

8E, < M{x)A'/ > Ey < M{y)A'/ > 

where A4 denotes any interpolating field for a bd meson. In eq. ( ^.4|) the time sep- 
arations —Xq and yo are to be taken large enough such that only the lowest state 
contributes, i.e. such that one obtains the on-shell matrix element. Analogously 
to the case of two-point functions, this condition is checked by looking for a joint 
plateau in xq and yo over some range of these variables. Such a plateau appears to be 
reached already at moderate —Xq, yo |7S|, |106| , PT| . Presumably, the reason is that the 



contributions of excited states to eq.( |5.4| ) are similar to the ground state contribution. 



Abada et al.[^ have computed the B-parameter for meson masses between 1.5 GeV 
and 3.5 GeV and with a cutoff a^^ ~ 3.7 GeV. Since the mass dependence is very weak, 
an extrapolation to rriB seemed justified. Also the B-parameter for the D meson was 
quoted: 

B^'^' = 1.05(8) (5.5) 

5|^^ = 1.16(7). (5.6) 
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The B-parameters turned out to depend only weakly on the mass of the light quark 
in the meson. For the ratio of the B-parameter of the strange meson to that of the 
nonstrange meson ref. |81| obtained (independently of the mass of the heavy quark) 



BpjBp^ = 1.02(2) . (5.7) 

We point out that the above quoted errors do not contain the O(a^) uncertainty in 
the renormalization (which includes mixing) of the operator. An estimate of the 0{a) 
terms in this quantity has not been given either. As discussed in the previous section, 
the latter uncertainty needs to be studied through systematic computations at different 
values of the lattice spacing. Before such a study is completed, the numbers quoted 
in this section have to be considered as rough estimates. 
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6. The Size of l/rrih Corrections to the Heavy Quark Limit 



An important nonperturbative question in heavy quark physics is to determine for 
which values of the heavy quark mass the predictions that one obtains in the heavy 
quark hmit, i.e. for mh — ^ oo, become accurate. Of course, this is not a question with 
a precise answer: if we have a quantity hke / with an expansion in the inverse heavy 
quark mass mj^^ 

f{ml') = t^\l + hml' + ...) , (6.1) 

the mass scale /i, that characterizes the size of the corrections is nonuniversal; it 
depends on the quantity considered. Nevertheless, it is of interest for applications of 
the HQET to see, whether - in examples - scales like /i are of the order of the size of 
the QCD scale Aqcd as it is frequently assumed. 

Fig. |12| gives a rough estimate /i ~ IGeV. Although a number of this order has been 
quoted in almost all recent publications [0, ^, ^, ^ , it is evident from our discussion 



in section ^ that /i has very large uncertainties. The main uncertainty is due to the 
present error of the static value. 

Hence, it is very useful to consider a quantity for which the heavy quark limit is known. 
Such a quantity was first discussed in refs. | |105| , |ST]]. Besides the pseudoscalar decay 



constant, one also considers the vector meson decay constant fy, which is convention- 
ally defined as 

< OKmV >= t^Fy^Ml/^/V2 . (6.2) 

Here, V^ denotes the (renormalized) heavy light vector current, |V^ > is a momentum 
zero vector meson state with polarization vector e^ and My is the mass of the vector 
meson. Due to the spin symmetry, the combination U{M) = FyFp/M, with M = 
{Mp + 3My)/4, becomes one in the heavy quark limit p3| : 



f/(M) = FyFp/M = (1 + ^a{M) + ...)(1 + Ui/M + ...) . (6.3) 

ovr 

Since U{M) is the ratio of two very similar quantities, the finite lattice spacing effects 
in U{M) may be reduced compared to the ones in / and also excited state contributions 
are expected to be less of a problem. Using masses around the mass of the D-meson, 
Baxter et al.[^ determined U{M) with the improved action. Interpreting their results 



in terms of f/i, one has f/i ~ 0.5 GeV in agreement with the rough results ||105| , ^ 
that were obtained with the Wilson action. 

So scales of about 0.5 GeV to 1 GeV are found for the coefficients of the l/rrih terms 
in the two known examples. This suggests that the heavy quark limit only gives 
the correct picture for charm quarks on a qualitative level, while it appears to be a 
quite reasonable approximation for 6-physics. In the important application of 6 ^ c 
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transitions, corrections of order 0(l/?7ic) may be of practical relevance. Although these 
conclusions are obtained entirely in the quenched approximation, it is to be expected 
that this qualitative result also holds true for full QCD. 

Since the l/rric-corrections are nonuniversal, the simple examples considered above 
give us only a rough idea about the quality of the approximations for different pro- 
cesses. In the end, the relevant quantities have to be determined directly from QCD. 



7. The B — B Threshold, String Breaking and Hybrid Mesons 



Following ref. |^3|, we now review, how to obtain information about the breaking of 
the QCD-string from the binding energy Mp of a static B-meson computed in the 
quenched approximation. Moreover, it is of interest to compare the B — B threshold 



with the lowest level of a hybrid h — h meson ||107| , |108|] in order to estimate whether 
such h — h mesons are broad or narrow resonances. This is done at the end of the 
section. 

In full QCD simulations, the breaking of the QCD string, i.e. the flattening of the 
heavy quark potential at large distances, has been sought after for some time. No 



effect was found in the most serious effort |58|. In the following, we estimate the 
distance r(,, where the full QCD potential flattens off using only quantities calculated 
in the quenched approximation |^ . 



Consider a large Wilson loop Vr(r, xq), with xq » r, in full QCD. It has a represen- 
tation in terms of the eigenvalues of the QCD-hamiltonian (or transfer matrix): 

W^(f,Xo) = Elcr(^1l'exp(-K(r>o) • (7.1) 

n>0 

Here, exp{—aVn{r)) are the eigenvalues of the transfer matrix in the corresponding 
charged sector of the Hilbert space: the states in this sector transform according to the 
3-representation at position and according to the 3-representation at position r under 
gauge transformations. The same states contribute in the spectral decomposition of 
the correlation function of static-light meson fields A4'^{x,xo) 

H{r,xo) = <M^{0,xo)[M-^{f,xo)]^M^{f,0)[M^{0,0)]^> 

= j:\^nirl\'eM-Vnir-)xo) . (7.2) 

ra>0 

The ground state potential V{r) = Vo(r) can therefore be called a static quark potential 
or a static meson potential. Physically, the first interpretation is sensible at relatively 
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short distances, whereas the second one is the appropriate language for large distances 
(compared to the confinement scale). This can even be put into a quantitative relation: 
we expect |c^(?^P << |cj5^(^P at small distances and |c^(^)P >> |cjf (r)p at large 
distances. 

Furthermore, at large distances, the potential will approach a constant up to nonlead- 
ing terms (Yukawa- type interactions), because the correlation function factorizes: 

log[H{r, xo)] = 2 \og[C'^'{xo)] + O(exp[-m,r]xo/r) . (7.3) 



Simulation results of full QCD [^] indicate that the QCD-potential is approximated 



rather well by the quenched potential out to relatively large distances, up to about 
r ~ 0.7 fm. At very large distances, on the other hand, we expect that H{r,t) is 
represented with some accuracy by the quenched approximation. We must switch 
from one correlation function to the other when using the quenched approximation, 
since we have to put the breaking of the string in by hand. Obviously, the quenched 
approximation does not give us much information about the intermediate regime. 

The (in r = \r\) asymptotic behavior of the correlation function H{f,t) is given by 
the mass in the static approximation. So one defines |^3| the string breaking distance 
n by 

Vin) = 2Mp, (7.4) 

with V{r) being the quenched potential including the self-energy term that cancels 
in eq.( [7.4| ). r^ defined in this way, gives an upper bound to the distance where the 



potential starts deviating significantly from the form of eg. (2.34) 



From the explicit parametrization of the potential eq.( |2.34| ) one can calculate r^ in 
lattice units: 




R, = (Mp - ^K.)/S + ^/ [{Mp - ^Vo)/n^ + ^ . (7.5) 



Extending the numerical evaluation of ref. [43|, we have included the data of refs. [35, 



and extrapolated rbvS linearly to the continuum (see fig. ITl). The resulting value 



in the continuum limit is (for vanishing light-quark mass) 

rbv^=3.5(2)(2), (7.6) 

where the first error is statistical and the second is an estimate of the uncertainty 
due to the extrapolation form (In this case, the a-dependence is rather strong). With 
^/o^ ~ 420 MeV, this corresponds to ri, = 1-7(2) fm. 

Such distances are difficult to reach in a calculation of the potential including dynam- 
ical fermions. It should be noted, however, that as a result of the investigation of 
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Figure 13: Continuum extrapolation of ri, for light-quark mass zero. Eq. (7.5) is used with data 
for Mp from refs. H H ||. 



the screening of the potential is expected to appear at distances that are only weakly 
dependant on the dynamical quark mass. The mechanism can hence be studied with 
relatively large (dynamical) quark masses. 

In the above presentation, we have translated the information that is present in the 
threshold energy 2Mp into the string breaking distance r;,, where the divergent self- 
energy of the static quark is cancelled. In a certain approximation, to be discussed 
below, the difference between the B — B threshold and the mass of a hybrid h — h 
meson [ll08| is also free of this divergence. As the position of the h — h energy level 
relative to the threshold determines the width of the h — h resonance, this quantity is 
an interesting observable. 

A hybrid h—b meson is a h—h bound state, with quantum numbers that do not occur in a 
naive potential model. Such "exotic" bound states are an interesting possibility due to 
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QCD. Describing the b-quarks through a nonrelativistic potential model, exotic states 
are possible through nontrivial angular momentum configurations of the gluon field 
which generates the potential. Such excited potentials have been computed from lattice 



gauge theory 1^0^]. The nonrelativistic Schrodinger equation with those potentials can 
be solved and thus one gets an estimate for the energy levels of hybrid 6 — 6 mesons. As 
the (excited) potential between static quarks is used, the energy levels contain (twice) 
the self energy of a static quark. The difference of an energy level to 2Mp is free of 
this divergent contribution and represents a well defined estimate of whether the state 
is above or below the open b threshold. Of course, the estimate is obtained within the 
quenched and the potential model approximations. 



Using the result of ref. ||107|| one estimates for the lowest exotic level in the b — b system 

[^i-+-n^o)]/v^~2.0 (7.7) 

and for the threshold 



[Ebb - V{ro)]/V^ = 2Mp - V{ro)/v^ = 2.5(4) 



(7.8) 



where we have arbitrarily normalized to V{rQ). So, the exotic state appears to be 
~ 200 MeV below the threshold. Adding the information that the wave function at 
the origin is very small due to the fiat potential ||107|| , the state is expected to have a 
very small width. Note, however, that even within the model used here, the uncertainty 
is so large that a coincidence of the threshold and the hybrid level is just at the edge 
of the error bar. 



Previously, the position of the BB threshold relative to the mass of the hybrid meson 
has been estimated in the following way ||107|| : The difference of Ei-+ and the mass 
of the lowest b — b bound state mr was estimated with the above model. Then the 
experimental value for 2mB — mr gives the position of the threshold relative to Ei-+ . 
In this way one arrives at the conclusion that the l""*" state is ~ 250 MeV above the 
threshold. 



Remaining consistently in the quenched approximation, the opposite seems to emerge, 
giving new hope for the existence of a narrow, interesting state. Further progress 
about this problem can most likely be made within the framework of non-relativistic 
QCD. 
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8. The Beauty Spectrum 



The spectrum of bound states of a b-quark with hght quarks represents one of the 
possible predictions of QCD. It has been explored to some extent using the static 
approximation. Estimates of the corrections of order 0{mi^ /itlb) do not exist and are 
in fact not of primary interest at this point, since the computations within the static 
approximation need considerable improvement. For any state, the self energy of the 
static quark must be cancelled to obtain a finite result. This is done by considering 
energy differences to the mass of the B-meson . 

Due to the spin symmetry in the static approximation 0, there is a degeneracy of 
vector- with pseudoscalar mesons and axialvector- with scalar mesons. We denote 
the remaining splitting between scalar and pseudoscalar by A5. Of further interest is 
the splitting between the A;, and the 5, denoted by Aa. Also the first radially excited 



state with splitting A2S of the B has been investigated! 95] and an estimate of the 



difference between the Bg and the B = B^ was given |BB]. We now comment on some 
technical points of these computations. 

The computation of the mass splittings between states with different quantum num- 
bers is straight forward in principle, i) One chooses (in addition to eq. ( |2.29| ) ) an 



interpolating field with the desired quantum numbers, ii) One searches for the ex- 
ponential decay of the correlation functions and computes the mass splitting as the 
difference, or one takes the ratio of the two correlation functions directly and searches 
for a plateau in the effective mass of the ratio. 

The scalar meson is a p-wave in the non-relativistic quark model. It is therefore natural 
to insert a p-wave smearing function into eq. ( p.28| ) in order to obtain an interpolating 
field for a scalar meson. For gauge covariant wave functions, this can be achieved 
by applying a covariant derivative to the symmetric wave function. This approach 



was tried in ref. E3], but the correlation function was found to be very noisy. The 



results quoted in ref. [|^ were obtained by replacing 7075 — > 1 in eq. (|2.29|) and using 



a spatially symmetric wave function, that means the lower components of the light- 
quark field are used to obtain the parity change relative to the pseudoscalar. Ref. [^ 
does obtain results with a p-wave smearing function in Coulomb gauge. Presumably, 
the reason is that in a fixed gauge, the detailed form of the smearing function can be 



chosen; ref. [^ uses the wave function of a semi-relativistic potential model for that 
purpose. 

The interpolating field for the A^ has been chosen as ||109| , ^ 

B'/{x)= Y: e,,,{h\x))l{u\x%{C^,)p,{d\x)r^ . (8.1) 

Here a,b,c (a, P, 7) denote color (Dirac) indices. The field B^^°'^{x) is to be interpreted 
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as an extreme di-quark trial wave function for the baryon: the two hght-quark fields 



are taken at the same point. Otherwise, B^ 



i,j. 



X 



amounts to a more general wave 



function ansatz. In ref. ||109|| , the di-quark option was chosen, while ref. [Q explored 



both possibilities. An interesting feature of the correlation functions of B'^'^{x) is that 
the di-quark interpolating field ;B^''°^(a;) is much less effective in exciting a A^ state 



than the field B^' {x), where all quarks have an independent spatial wave function 1 43 



Ref. AA[GeV] As[GeV] AsJGeV] M b - M^jMeV] a-^[GeV] 

ra o.72t^:in.i3 



d5 



m 



' J--13-16 



0.6 



0.35 
^0.4 



0.4 



70 - 140 



1.9 

1.3-2.7 

1.6 

1.9 



Table 4: A compilation of mass splittings as they were obtained in the static approximation. The 
last column gives the value of the cutoff. 



Determining the energy of a (radially) excited state is a more difficult task. For that 
purpose, a matrix correlation function needs to be considered. From such a matrix 
correlation function, a quantity which is analogous to the local mass (cf. sect. |2.5|) can 
be constructed and it was shown that this quantity converges exponentially (in xq) to 
the mass splitting [ [L10|| . Thus the situation is - at least in principle - the same as for 
the splitting between two states with different quantum numbers. 



The computation of the 2S - IS splitting in ref. [95] has not exactly been done in 
this way. Rather, starting from a 2 x 2 matrix correlation function (where the two 
smearing functions are again obtained from the semi-relativistic potential model and 
should lead to relatively good interpolating fields for the lowest two levels), an optimal 
wave function was found for the ground state first. Then the correlation function was 
projected onto the part orthogonal to that approximate ground state. The excited state 
mass was determined from that projected correlation function. Since there always 
remains a small contribution from the exact ground state in the projected correlation 
function, one will at very large Xq determine again the mass of the ground state. 

we do not think 
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Although this is in principle a problem of the computation ref. 

that this effect is numerically important at the moment. Furthermore, it can easily be 

corrected in the future. 



The main feature of all the computations [|T0^, ^, ^ of the mass splittings is that the 
plateaus in the effective masses are less convincing than for the pseudoscalar state. So 
the results are less precise and it has not been possible to study systematic errors due 
to a-effectsQ. Therefore, these estimates are qualitative at the moment. Where this 
was done in the literature, we will quote numbers with error bars below, but we stress 



^^ Finite volume effects have been studied within a model for As and Ao 
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that they do not include a reahstic estimate of systematic errors. The results of the 
various investigations are listed in table §. 

In addition to the observables that we considered so far, Bochicchio et al.[ |109|| esti- 
mated the vector - pseudoscalar splitting at order l/nih. The computation is done 
starting from the static approximation and including the aB term of eq. (|2.36|) per- 



turbatively. The result at a"! = 1.9 GeV is ml - m], = (0.19l|j:[}^l|j.o7)GeV^ which 
is significantly below the experimental iv?g, — iv?g ~ m?]j, — m]-) ~ 0.55 GeV^. Lack- 
ing estimates of the systematic errors, it is premature to speculate from where this 
difference originates. It poses an interesting problem for future investigations. 
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9. Further Lattice Investigations 



There are a number of interesting investigations of other physical observables in the 
context of HL hadrons. Let us mention just the ones which we beheve have the poten- 
tial to contribute to the analysis of experimental data as mentioned in the introduction. 



Semileptonic decay form factors for /^-decays have been studied [|T2|, |1^ and ex- 
ploratory studies of i?-decays have been performed [p!3|. Computations of the Isgur- 



Wise-function are attempted! 1 11 1 and the nonperturbative amplitude for the decay 



B -^ K*^ has recently been estimated [pT 

We mention these investigations for completeness. Compared to the leptonic decay 
constants there is scarcely a study of the systematic uncertainties available at the 
moment. In particular, the size of finite lattice spacing effects are not known at 
present, but we would like to point out that the feasibility of such computations is 
clearly demonstrated. Thus, there is considerable potential in these computations. An 
initial understanding of the systematic uncertainties of these quantities should develop 
within the next year or two. After that stage, these computations might represent 
predictions of a well-tested model, quenched QCD. As such, they might already be a 
valuable contribution to the analysis of experimental data. 

The necessary computations in full QCD will be facilitated by what can be learned from 
the quenched approximation. Nevertheless, it is not possible to predict exactly when 
first principle calculations of these hadronic matrix elements, i.e. reliable calculations 
in full lattice QCD, can be done. 
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10. Summary 



We have given an overview of the present status of QCD lattice simulations that 
involve the b quark. The discussion was centered around the most easy quantity, 
the vacuum-to-one- meson matrix element. That matrix element gives us the leptonic 
decay constant and in the case of the i?-meson it is the primary unknown in the 
process of extracting the CKM angles from experimental data. 

Unfortunately, the i?-meson is particularly difficult to treat in a lattice simulation. 
This is simply due to the large gap in the relevant scales in the problem. These scales 
are i) the confinement scale which determines the physical size of a 5-meson and 
thus is the scale that is relevant for finite size effects and ii) the mass of the 6-quark 
itself. The latter needs to be small compared to the cutoff such that the quark can 
propagate without large distortions due to the finite cutoff. These two scales cannot 
be accommodated on today's lattices. 

A way out is to consider the observables of interest in two unphysical regimes, obtaining 
the physical observables through a matching of the two. The first regime is the static 
limit rrib ^ oo. In this limit, the scale rrih itself becomes irrelevant and results can 
- in principle - easily be obtained. The second regime is ruh ~ rric, where one is 
at the edge of the possibility of treating the heavy quark correctly while still keeping 
finite size effects small. We have emphasized that it is very important to systematically 
perform simulations with different values of the cutoff a~^ and fixed value of rUh in this 
regime and extrapolate the results to the continuum limit a — ^ 0. The uncertainties 
due to the finite cutoff can only be taken into account in this way. The necessity to 
extrapolate a — > introduces quite large errors, even in the quenched approximation. 
Nevertheless, it is straightforward to reduce these uncertainties using the knowledge 
which has been acquired over the last years. 

The precise studies are still restricted to the quenched approximation. It is quite inter- 
esting, however, that a systematic comparison of full QCD results with the quenched 
approximation in the range rUh ~ rric shows no effect of sea quarks with a mass of 
about TTis- 

Other quantities like the B-parameter and mass splittings have not been studied very 
systematically yet. Particularly concerning the B-parameter, this deficiency needs to 
be filled. 

In addition, we explained how one can get information on the breaking of the QCD- 
string from simulations in the quenched approximation by comparing the open B 
threshold in the static approximation with the static potential. 



Renormalization of Vector Currents 



For HL-mesons the axial vector current and the vector current have been normahzed 
in different waysPB], Bl, P3|, as discussed for the axial vector current in sect. P]^. In 
the case of the local vector current, the relativistic normalization is 



VfJ\ 



Zvigl.Kf.Kf,) qf{x)-f^qf 



with the one-loop expression ||3T 



Zvigl Kf, Kf>) = 1 - 0.17408 g^ + 0{g^) + 0(amf) + 0{amf) . 
The non-relativistic normalization of Kronfeld, Lepage and Mackenzie reads 



V^^f'ix) 



~Zv{glKf,Kf 



2K, 



X] = Z:v{Qn.J^f.J^f'} Qf[X]1„Qf'[X] (Al) 

(A2) 

- 3^^ J^^ - 3m qf{x)-i^qf>{x) , (A3) 

Zvig'o, Kf, Kf,) = 1 - 0.0656 f + 0{g^) + 0{amf) + 0{amf,) . (A4) 

It is claimed in the literature that in the region where arrif is of order one, the nor- 
malization eqs. (X3,XiD (and eqs. ( 2.22 . 2!25| )) represent the "correct" normalization 
of the currents. Here, "correct" has to be interpreted as being the normalization with 
small lattice artifacts. 

The partially conserved vector current 

V^'f'{x) = -[g7(x)(7^ - l)f/^(x)g/(x + afi) + qf{x + afi){l + lp,)Ul{x)qf{x)] (A5) 



with 



has the same normalization in the two approaches, however [|1 13 



We consider [|I3|, |113|| the ratio 



RviKf,K 



f 



J:,<qf'{x)-fkqf{x) V/'^'(0)> 

Ex < Qf'{x)-fkqf{x) g/(0)7feg/'(0) > 



(A6) 



Away from the short distance lattice artifacts this gives one definition of the renor- 
malization constant of the local current. In the relativistic normalization one obtains 
for this ratio 

Rv{Kf,Kf,) = Zv{go,Kf,Kf>) = l-O.UAOSf + 0{f) + 0{amf)+0{amf,) , (A7) 

while the non-relativistic normalization predicts a strong mass dependence 



Rv{Kf,Kf,) = (1-0.0656^ 



s2\ 



2Kf 



— 3m 



2Kf> 



3m + 0{g'^) + 0{amf) + 0{amf 

(A8 
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In fig. n, the Monte Carlo data for Ry is compared to tlie above two expressions. It is 
seen that at small m = l/{2Kf) + l/{2Kf>) — 1/ K^ the predictions are essentially the 
same. At growing quark masses, the prediction of the non-relativistic normalization 
deviates much further from the Monte Carlo results than the one from the relativistic 
normalization. In fact, the latter describes the mass-dependence very well. One also 
sees that the difference between the 1-loop expression for Ry and the nonperturbative 
results decreases significantly as the lattice spacing is reduced (from the right figure 
to the left, the lattice spacing changes by a factor ~ 0.6). 





Figure 14: Ry compared to eq. ( [A7| ) (dashed line) and eq. (^ 
is for /3 = 6.4 111 and the right hand figure for /3 = 6.0 |ll|. 



1) (dotted hne). The left hand figure 



Finally, let us mention that for the improved action, one obtains a much better agree- 
ment of Ry with 1-loop perturbation theory. This and the fact that other nonper- 
turbative definitions of Zy give significantly different results pO|, ^ , suggests that the 
difference between the dashed curve and the data in fig. TA is dominantly an 0{a) 
lattice artifact rather than an effect of truncating perturbation theory at that order. 
Indeed, the numbers are in agreement with an almost linear correction in a, but does, 
of course, not rule out more complicated corrections. 



This example shows that the nonrelativistic normalization, in general, does not re- 
duce lattice artifacts that are due to large quark masses mja ~ 1. Rather, in the 
example, the lattice artifacts are increased compared to the relativistic normalization. 
We mainly conclude from this example that lattice artifacts have to be removed by 
systematically performing the limit a — > 0. 
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